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Abstract 

We review the status of bottom quark physics at the CDF experiment. 
The measurements reported are based on about 110 pb _1 of data collected at 
the Fermilab Tevatron pp Collider operating at y/s = 1.8 TeV. In particular, 
we review results on B hadron lifetimes, measurements of the time dependence 
of B°B° oscillations, and a search for CP violation in B° — > J/ipKg decays. 
Prospects for future B physics at CDF in the next run of the Tevatron Collider 
starting in the year 2000 are also given. 
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We review the status of bottom quark physics at the CDF experiment. The measure- 
ments reported are based on about 110 pb — 1 of data collected at the Fermilab Tevatron 
pp Collider operating at y/s = 1.8 TeV. In particular, we review results on B hadron 
lifetimes, measurements of the time dependence of B°B° oscillations, and a search for 
CP violation in B° — » J/tpK^ decays. Prospects for future B physics at CDF in the 
next run of the Tevatron Collider starting in the year 2000 are also given. 

1. Introduction 

In 1977, the bottom quark was discovered as a resonance in the dimuon invariant 
mass spectrum in 400 GeV proton-nucleus collisions at Fermilato. Soon after the 
discovery of this new bb bound state with a mass of about 9.5 GeV/c 2 , the T res- 
onances were confirmed in e + e~ collisions at the DORIS storage ring at DESYi. 
Today, more than 20 years later, all lowest mass bound states containing a b quark 
have been discovered and are experimentally well established. The pseudoscalar 
B meson states a are B° = \ bd), B + = \ bu), B% = \ bs), and B+ = \bc), while 
the A° — | bdiL) is the b baryon ground state with lowest mass. The rest masses of 
the B hadronsu are between 5.3 GeV/c 2 and 6.4 GeV/c 2 , approximately six times 
the mass of a proton. The lowest lying B hadrons decay via the weak interaction. 

B hadrons play a special role among hadrons. The heaviest quark, the top 
quark, decays weakly into a real W boson and a b quark before it is able to form a 
meson with another antiquark through the strong interaction. Therefore, hadrons 
containing a b quark are the heaviest hadrons experimentally accessible. The prin- 
cipal interest in studying B hadrons in the context of the Standard ModelS arises 
from the fact that B hadron decays provide valuable information on the weak mix- 
ing matrix, the Cabibbo-Kobayashi-Maskawa (CKM) matrixd!. In fact, B decays 
measure five of the nine CKM matrix elements: V c b, Vub, Vtd, Vts, and Vj&. The 
future interest in B physics certainly lies in the study of CP violation in the system 
of B mesons. 

Traditionally, B physics has been the domain of e + e~ machines operating on 
the T(45) resonance or the Z° pole. But the UA 1 collaboration has already shown 

a Throughout the paper, unless otherwise noted, references to a specific charge state are meant to 
imply the charge-conjugate state as well. 
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that B physics is feasible in a hadron collider environment (for a review see Ref .0) . 
Although B physics did not play a significant role in the considerations for the 
original 1981 technical design report of the Collider Detector at Fermilab (CDF), 
several features in the CDF design are advantageous for the studies of B decays. 
These features include a large solenoidal magnetic tracking volume, a well seg- 
mented calorimeter for the detection of electrons, and muon chambers that allow 
low momentum muon detection. However, the device that made B physics possible 
at CDF in a competitive way, allowing for a broad B physics program, is a silicon 
micro- vertex detector, installed in 1992. 

There are several motivations for pursuing B physics at the Fermilab Tevatron 
pp Collider operating at i/s =1.8 TeV. The primary reason is demonstrated by 
comparing the B production cross section in e + e~ collisions, which is about 1 nb 
at the T(45) and about 6 nb at the Z° pole, to the large b quark production cross 
section at a hadron collider. At the Tevatron, ab is ~ 50 /A> within the central 
detector region of rapidity less than one. This is a huge cross section which results 
in about 5 • 10 9 bb pairs being produced in 100 pb _1 of data. However, the total 
inelastic cross section at the Tevatron is still about three orders of magnitude larger 
than the b cross section. This puts certain requirements on the trigger system used 
to find B decay products, as will be further discussed later. 

In addition to high rates of B hadrons, the Tevatron offers other features worth 
noting. First, in contrast to an e + e~ machine operating at the T(4S'), all B hadron 
species are produced at a hadron collider. Second, the transverse momentum spec- 
trum for B hadrons scales with the B mass and is significantly harder for heavy 
hadrons than that for light hadrons. As a result, the B hadrons are Lorentz-boosted 
at all rapidities, including the central detector region where the production rate is 
the highest with an average transverse B momentum around 4-5 GeV/c. Third, 
the hard B hadron momentum spectrum can be exploited to improve the signal 
to background ratio in finding B decay products. For all momenta, b production 
accounts for about 0.2% of the total pp inelastic cross section, while at high mo- 
menta the ratio of b jet to inclusive jet production is close to 2%. Finally, the initial 
pp state is a CP symmetric state where we expect equal rates of B and B hadrons 
to be produced, at least in the central detector region. 

In this article, we review recent B physics results at CDF concentrating on B life- 
time measurements, proper time dependent measurements of B°B° oscillations, and 
the search for CP violation in B° — > J/ipKg. The measurements reported here 
are based on about 110 pb _1 of data collected at the Fermilab Tevatron Collider. 
Prospects for future B physics at CDF in the next run of the Tevatron starting in 
the year 2000 are also discussed. The outline of this paper is as follows: In Section 2, 
we give a brief overview of heavy quark production in pp collisions. In Section 3, 
we describe the experimental environment including the Tevatron Collider, as well 
as the CDF detector. We focus on the collection of B physics datasets at CDF 
in Sec. 4, emphasizing the CDF trigger scheme. In Section 5, we highlight several 
features of B physics in a hadron collider environment and describe some of the 
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ways B decays are studied at CDF. Section 6 is devoted to the measurements of 
B hadron lifetimes, where the lifetimes of all weakly decaying B mesons as well as 
the A° baryon are measured at CDF. We then review several proper time depen- 
dent measurements of B°B° oscillations in Section 7 and discuss various B tagging 
methods to identify the B flavour in hadronic collisions. The search for CP viola- 
tion using the current data set of B° — > J/ipKg decays is summarized in Section 8. 
An outlook for future B physics at CDF in the next run of the Tevatron Collider, 
starting in the year 2000, is given in Sec. 9. Finally, we offer our conclusions in 
Section 10. 

2. Heavy Quark Production in Hadronic Collisions 

In this section, we give a short introduction to heavy quark production in pp col- 
lisions. With the expression "heavy quark" we will mainly refer to b quarks. The 
discussion will be rather qualitative and cannot serve as a complete review of heavy 
quark production in hadronic collisions. It is meant to give an idea of the main 
issues, introducing some of the nomenclature ofterLjjsed in the literature. The in- 
terested reader is referred to reviews in referencesKM. Here, we first discuss the 
parton model and parton distribution functions. Then, we summarize heavy quark 
production in lowest order QCD, as well as next-to-leading order QCD calculations 
and discuss hadronization of heavy quarks. Finally, we briefly compare theoretical 
predictions of b quark production to measurements at CDF. 

2.1. Parton model and parton distribution functions 

In a static picture, a proton is a bound state of three quarks | uud) with a radius 
of about 1 fm. However, in a hadronic collision such as that found at the Teva- 
tron, a proton is better characterized as a beam of free partons: Three constituent 
quarks (valence quarks), virtual gluons, and quark-antiquark pairs (sea quarks )n. 
The different partons don't necessarily divide up the beam energy equally. The 
distribution of partons within the proton is described by the so-called parton dis- 
tribution function (PDF) F°"{x, Q 2 ) which is the number density of parton i (quark 
or gluon) carrying the momentum fraction x of the hadron a (proton or antiproton) 
when probed at momentum transfer Q 2 . 

Information on parton distributions comes from measurements of deep inelastic 
lepton or nucleon scattering, such as ep, eN ', /xiV, or vN. Once parton distributions 
have been measured at some value of Q 2 (normally at low Q 2 ), and the running cou- 
pling constant a s (Q 2 ) of the strong interaction has been determined, QCD permits 
to compute the parton distributions at high er values of Q 2 based on the formalism 
provided by the Altarelli-Parisi equationsEil and the evolution calculations pioneered 
by Dokshitzer, Gribov, and LipatovtL3. Modern representations and predictions of 
parton distribution functions, used for comparison with hadron collider data, come, 
for example, from the Martin-Roberts-Stirling (MRS) grouptSJ or frorn a collab- 
oration of theorists and experimentalists called the CTEQ collaboration^. 
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Fig. 1. (a) Q 2 evolution of the momentum-fractions carried by the various parton species in the 
proton for the CTEQ4 parton distributipiiscj ■ (b) Total production cross sections for charm, bot- 
tom, and top quark pairs in pp collisions! 10 ! as a function of the hadronic centre-of-mass energy y/s. 



As an example, Figure |l|a) illustrates the flavour content of the proton as mea- 
sured by the momentum fraction L dxxFi(x, Q 2 A, carried by each parton species i, 
as obtained for the CTEQ4 parton distributional^. Gluons carry about half of the 
proton's momentum almost independently of Q 2 , while the momentum is shared 
more and more equally among the quark and antiquark flavours as Q 2 increases. 

2.2. QCD predictions for hard scattering in pp collisions 

In a pp interaction at the Tevatron, heavy quarks Q are produced in the hard 
collision of two partons, one from each hadron. In terms of the rapidity y, defined as 



2 \E- Pz 

and transverse momentum p t , the relativistically invariant phase space volume ele- 
ment of the final state heavy quark with four- momentum (E,p x ,p y ,p z ) is 

~dydV (2) 

The differential partonic cross section &ij per invariant phase space volume for the 
production of a heavy quark in a given parton-parton subprocess can be written as 

d 3 ^ Ed 3 &ij 

d^p/E = d 3 p ^ XiPa ' XjPb ' P]mQ ' A>fJ ' R ' fJ ' F ^ 

where Xip a and XjPb are the momenta of the incoming partons i and j, p is the 
momentum of the outgoing heavy quark, while p a and pb are the momenta of the 
colliding hadrons. The remaining variables are parameters of the theory: tuq is 
the mass of the heavy quark, A determines the coupling strength of the strong 
interaction expressed through a s , (1r is the renormalization scale, which is related 
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to the energy scale used in the evaluation of a s , and /zp is the factorization scale 
used in the evolution of the parton densities. In QCD, the partonic cross section can 
be expressed as a perturbative expansion in a s , provided mg is sufficiently large. 

The cross section for heavy quark production in the collision of two hadrons 
(a + b — > Q + X) is then obtained by convoluting the partonic cross section above 
with parton distribution functions Fi(x,Q 2 ), where it is common to equate the 
momentum transfer Q 2 = m 2 g + p\ with /i 2 = [i 2 R = [i 2 F . Note, the uncertainty 
on the choice of the /i-scale, which is not extractable from the data, is one of 
the large sources of uncertainties in QCD predictions of heavy quark production. 
The perturbative QCD formula for the inclusive production of a heavy quark in a 
hadron-hadron collision is then given as 

—^T~ = Y1 J dxzdxj f ^3°^ (xiP a , XjPb,p; m Q , A, \ F°'(x ll ^ 2 )F b j {x j ^ 2 ). 

ij 

(4) 

The corrections to Eq. (Q) are suppressed by powers of the heavy quark mass. 
Finally, the total cross section for the production of a heavy quark is obtained 
by integrating Eq. (Q) over momentum p. As a further illustration, Figure [l]b) 
compares the total production cross sections for charm, bottom, and top quark 
pairs in pp collisionst2l as a function of the hadronic centre-of-mass energy y/s. 
Note the different units used for top quarks compared to charm and bottom. 

2.2.1. bb production in leading order QCD 

The leading order (LO) a 2 diagrams for heavy flavour production are the 2^2 
processes of gluon-gluon fusion g + g — > Q + Q as shown in Fig. ||a-c) and quark- 
antiquark annihilation q + g — ► Q + Q displayed in Fig. |^d) , respectively. In the 
latter process, the QQ pair is always in a colour octet state while in g g — > Q + Q 
both colour singlet and octet are allowed. The gluon-gluon fusion process is the 
dominant production mechanism for b quarks at the Tevatron, while top quarks 
are mainly produced from the quark-antiquark annihilation process. The lowest 
order 0(a 2 ) matrix elements in the a s expansion and the cxoss sections for these 
processes have been available in the literature for some timal3. 

We briefly summarize some of the phenomenological consequences of the 0{a 2 ) 
processes. First, we note that the partonic cross section <x is proportional to a 2 /m 2 Q. 
The average transverse momentum of the heavy quark grows approximately with 
its mass (pt(Q)) ~ m Q, meaning that the average b quark transverse momentum 
is about 4-5 GeV/c. In addition, the pt distribution falls rapidly to zero as pt 
becomes larger than the heavy quark mass. Furthermore, the heavy quark and 
antiquark are produced back-to-back in the parton-parton centre-of-mass rest frame 
and are correspondingly back-to-back in the plane transverse to the colliding hadron 
beams. In addition, the rapidity distribution of the QQ pair has a typical bell shape 
becoming wider and flatter as the partonic energy grows. This means, heavy quark 
production is larger in the central region falling at higher rapidities. The rapidity 
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Fig. 2. Lowest order 0(a^) diagrams for heavy quark production through (a)-(c) gluon-gluon 
fusion and (d) quark-antiquark annihilation. 

difference between the Q and Q tends to be of order unity. But for a fixed value of 
Pt, the production rate is highly suppressed when the rapidity difference becomes 
large, which means that Q and Q tend to be produced with similar rapidity. 

2.2.2. bb production in next-to-leading order QCD 

The next-to-leading order (NLO) terms in the a s expansion were originally consid- 
ered 'corrections' to the leading order terma^ But, it was soon recognized that 
the higher order corrections could be largetS It was noticed that the process 
.9 + 9 ~~ * .9 + 9 with g — > Q + Q, which is formally of order a s , can be as im- 
portant as the lowest order processes because the cross section for the production 
of gluons g + g — > g + g is about a hundred times larger than the LO cross section 
for the process g + g — ► Q + Q. The complete calculations of next-to-leading order 
corrections to hea vy q uark production in hadronic collisions have been performed 
by several authorst3E3. It has been shown that the O(a^) terms are actually larger 
than the lowest order processes for b and c quark production, if the hadron-hadron 
centre-of-mass energy is much larger than mg which is the case for the Tevatron. 

Examples of order diagrams are displayed in Figure |. These processes in- 
clude real emission matrix elements (Fig. ^a) and the interference of virtual matrix 
elements with the leading order diagrams shown in Fig. |j^b). Other NLO contribu- 
tions come from gluon splitting diagrams (Fig. ^Jc) or the flavour excitation process 
as displayed in Fig. ||d) . In the gluon splitting case, the probability to find a heavy 
quark from a gluon with large p% has a logarithmic increase. Phcnomcnologically, 
the QQ pair is produced close in phase space and will often appear as a single 
jet. In the flavour excitation process, the heavy quark is considered to be already 
present with a certain heavy quark density in the incoming hadron. It is excited 
by the exchange of a gluon with the other hadron and appears on mass-shell in the 
final state. In the case of a flavour excitation process, only one of the quarks from 
the QQ pair is usually at high p t . 

In summary, next-to-leading order processes are an important contribution to 
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Fig. 3. Examples of diagrams for heavy quark production at next-to-leading order: (a) Real 
emission diagrams, (b) virtual emission diagrams, (c) gluon splitting, and (d) flavour excitation. 

b quark production in pp collisions at =1.8 TeV. These processes are the source 
of bb pairs that can be close together in phase space and do not necessarily appear 
in the central detector region as compared to the LO processes which result in back- 
to-back quark pairs. However, there is evidence (see Sec. 2.4.) that NLO processes 
are not sufficient to obtain accurate estimates since large /i-scale dependences are 
still present. These are an indication that higher order corrections could be large. 

2.3. Hadronization of b quarks 

Once b quarks are produced through the initial hard scattering, the process of 
forming B hadrons follows and is called hadronization or fragmentation. It is a 
low Q 2 process and is beyond the reach of perturbative QCD calculations. The 
hadronization process is therefore described by semi-empirical models inspired by 
theory. A commonly used approach is the string fragmentation modeEI. In a 
naive picture, we can imagine a "cloud" of gluons acting as a string between the b 
and b quark pair. As the quark and antiquark separate, the string stretches until 
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it breaks and a new qq pair is created out of the vacuum to form the new ends 
of the string. These new strings also stretch and break, producing more quark- 
antiquark pairs until there is no longer sufficient energy available to generate new 
qq pairs. The particles produced in this hadronization process, along with the 
B hadron, are usually referred to as fragmentation particles. The fraction of the 
initial b quark momentum transfered to the B hadron is commonly described by a 
so-called fragmentation function, suggested, for example, by Peterson et alE3 for 
the case of a heavy quark Q forming a hadron together with a light quark q: 

dAr oci.fi-i-^y 2 . (5) 



dz z \ z 1 

Here, is the so-called Peterson fragmentation parameter, related to the ratio of the 
effective light and heavy quark masses £& ~ {jriq/mg) . The variable z is originally 
defined as z = (E + p\\)Qq/{E + pq) where p|| is the projection of the momentum 
of the hadron onto the direction of the heavy quark before hadronization. This 
variable is often approximated with experimentally better accessible parameters 
such as x p —p/pmax or xe — -Ehadron/£ : bcam used in e + e~ annihilation. 

2.4. Comparison of heavy quark production with CDF data 

To conclude our review of heavy quark production, we compare in Figure ^, taken 
from Ref.EI, various CDF measurements of b quark production with next-to-leading 
order predictions. The distribution most commonly studied by hadron collider ex- 
periments is the b quark differential p t spectrum, integrated above a given p t thresh- 
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CDF Integrated Luminosity -- Run I (pb-1) 

End Run 1 Feb. 19. 1996 




Fig. 5. Time profile of the luminosity delivered by the Tevatron Collider and accumulated by the 
CDF experiment throughout Run I. 

old p™ m and within a fixed rapidity range which is usually set to j/ ma x = 1 at CDF: 
a(p t >pf n )= I dy f dPtJ^T- (6) 

Different techniques, involving primarily CDF's silicon vertex detector to improve 
background rejection, are used to obtain the data points from, e.g. high-p t J/tp 
mesons reconstructed through fi + fi~ or high-p t leptons accompanied by a nearby 
charm meson (D°). The measurements are compared to next-to-leading order QCD 
calculations which agree very well in shape with the data. However, the data are 
higher by a factor of almost three than the default prediction based on /i = /iq = 
y/p% + m%. Using variations on, for example, m&, /i, and the parton distribution 
function, theory is able to accommodate the data as indicated by the upper curve 
in Fig. ^. Since the production o£J> quarks and B hadrons is not subject of this 
review, we refer to the literatura§cJ for a further discussion of this issue. 

3. Experimental Environment 

In this section, we summarize the experimental environment beginning with the 
Tevatron Collider. We then describe the CDF detector, emphasizing the device that 
made a broad B physics program possible at CDF, the silicon micro- vertex detector. 

3.1. Tevatron collider 

The Fermilab Tevatron, with a circumference of 6.28 km, is a proton-antiproton 
collider operating at a centre-of-mass energy of y/s = 1.8 TeV. The Tevatron deliv- 
ered its first physics collisions in 1987. In this article, however, we will concentrate 
on the data taking period referred to as Run I which started in May 1992 and ended 
in Feb. 1996. During that time, the Tevatron operated with six bunches of protons 
and six bunches of antiprotons crossing every 3.5 us at CDF's interaction region. 




Fig. 6. Schematic view of the CDF detector. 



Figure [5j shows a time profile of the delivered and accumulated luminosity at 
CDF throughout Run I. During this period, a total integrated luminosity of about 
f 75 pb^ 1 was delivered to the CDF experiment. During Run I, the highest instanta- 
neous luminosities reached, were around 2.5 • 10 31 cm s . The Runl running pe- 
riod was divided up into a Run la, from May 1992 through June 1993, a Run lb, from 
Dec. 1993 to July 1995, and a Runic, from Dec. 1995 to Feb. 1996. The collected 
integrated luminosities of data used for physics analyses were about 19.3 pb _1 for 
Run la, approximately 90 pb" 1 for Run lb, and about 10 pb _1 for Runic. However, 
the Run Ic data taking period was dedicated to the accumulation of very special- 
ized trigger datasets and is usually not included in Runl physics analyses, resulting 
in about 110 pb -1 of data used for physics results from Runl. All measurements 
presented in this paper refer to this Run I luminosity, unless otherwise noted. 



3.2. CDF detector 

The CDF detector is a multi-purpose apparatus designed to study 1.8 TeV pp col- 
lisions produced by the Fermilab Tevatron Collider. It has both azimuthal and 
forward-backward symmetry. A superconducting solenoid, 4.8 m in length and 
1.5 m in radius, generates a 1.4 T magnetic field containing tracking chambers used 
to detect charged particles and measure their momenta. Surrounding the solenoid, 
sampling calorimeters measure electromagnetic and hadronic energies of jets, elec- 
trons, and photons. Outside the calorimeters are drift chambers for muon detection. 

A schematic view of the CDF experiment can be found in Figure p while a side- 
view quarter cross section of the CDF detector is displayed in Fig. 7|. The CDF 
experiment uses a coordinate system with the z-axis along the proton beam direc- 
tion, the y-axis pointing vertically upwards, and the x-axis pointing horizontally 
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Fig. 7. Side-view quarter cross section of the CDF detector. For a full view, the image is to be 
reflected to the left and below. The pp interaction point is at the lower-right corner of the figure. 



out of the Tevatron ring. Throughout this article, (p is the azimuthal angle, 9 is the 
polar angle measured from the proton direction, and r is the radius perpendicular to 
the beam axis. The pseudorapidity r\ is defined as r\ = — ln[tan((9/2) ]. The trans- 
verse momentum p t is the component of the track momentum p transverse to the 
z-axis (pt = p ■ sin 9) while E t — E ■ sin 9 with E being the energy of the calorimeter 
cluster. A more complete description of the CDF detector can be found elsewheretfl. 
We summarize here only those detector features most relevant to B physics. 

Three devices inside the 1.4 T solenoidal magnetic field are used for the track- 
ing of charged particles: The silicon vertex detector (SVX), a set of vertex time 
projection chambers (VTX), and the central tracking chamber (CTC). The SVX is 
described in more detail in Section 3.2.1. The VTX, which is located outside the 
SVX up to a radius of 22 cm, reconstructs track segments in the r z-plane up to 
1 77 1 < 3.25. The VTX is used to determine the z-position of the primary interaction 
vertex with a resolution of about 0.2 cm on average. 

Surrounding the SVX and VTX is the CTC, located between radii of 30 cm and 
132 cm. The CTC is a 3.2 m long cylindrical drift chamber that contains 84 layers 
of sense wires grouped into nine alternating super-layers of axial and stereo wires 
with a stereo angle of 3° . The CTC provides three-dimensional tracking and covers 
the pseudorapidity interval \r]\ less than about 1.1. The outer 54 layers of the CTC 
are instrumented to record the specific ionization dE/dx of charged particles. 

Outside the solenoid are electromagnetic (CEM) and hadronic (CHA) calorime- 
ters (| 77 1 < 1.1) which employ a projective tower geometry back to the nominal 
interaction point with a segmentation of A77 x Aip ~ 0.1 x 15°. The sampling 
medium is composed of scintillators layered with lead absorbers in the electromag- 
netic calorimeter and steel in the CHA. The energy resolution for the CDF central 
calorimeter is a (E t )/E t = [(13.5%/ y/K) 2 + (2%) 2 ] 1 / 2 for electromagetic showers 
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and a{E t )/E t = [(50%/y^) 2 + (3%) 2 ] 1/2 for hadrons with E t measured in GeV. A 
layer of proportional chambers (CES), with wire and strip readout, is located six ra- 
diation lengths deep in the CEM calorimeters approximately near shower maximum 
for electromagnetic showers. The CES provides a measurement of electromagnetic 
shower profiles in both the ip- and z-directions. Proportional chambers located be- 
tween the solenoid and the CEM comprise the central preradiator detector (CPR) 
which samples the early development of electromagnetic showers in the material 
of the solenoid coil, providing information in r-ip only. Finally, plug and forward 
calorimeters instrument the region of 1.1 < \rj\ < 4.2. These consist of gas propor- 
tional chambers as active media, while lead and iron are used as absorber materials. 

The central calorimeters act as hadron absorber for the muon detection system. 
Four of its layers of planar drift chambers (CMU) are located beyond the central 
calorimeters. The CMU system covers \rj\ < 0.6 and can be reached by muons with 
Pt in excess of 1.4 GeV/c. To reduce the probability of misidentifying penetrating 
hadrons as muon candidates in the central detector region, four additional layers 
of drift chambers (CMP) were added in 1992 and are located behind 0.6 m of steel 
outside the CMU system. Approximately 84% of the solid angle for \r/\ < 0.6 is 
covered by the CMU detector, 63% by the CMP, and 53% by both. To reach 
these two detectors, particles produced at the primary interaction vertex, with a 
polar angle of 90°, must traverse material totaling 5.4 and 8.4 pion interaction 
lengths, respectively. An additional set of muon chambers (CMX) is located in the 
pseudorapidity interval 0.6 < |?7| < 1.0 to extend the polar acceptance of the muon 
system. Approximately 71% of the solid angle for 0.6 < | ^-7 1 < 1.0 is covered by the 
free-standing conical arches of the CMX. Finally, the forward muon system (FMU) 
is a magnetic spectrometer consisting of three planes of drift chambers surrounding 
two 1 m thick iron toroids located ±10 m from the interaction point. 

3.2.1. CDF silicon micro-vertex detectors 

Surrounding thp_1.9 cm radius beryllium bcampipe is the CDF silicon micro- vertex 
detector (SVX)@ originally installed at CDF in 1992. The SVX is 51 cm long and 
consists of two identical cylindrical modules which meet at z = 0, with a gap of 
2.15 cm between them. A sketch of one of these modules is shown in Fig. |8^,). Both 
SVX modules consist of four layers of silicon micro-strip detectors segmented into 
twelve 30° wedges. The basic detector element is called a ladder and is shown in 
Fig. ||b). There are 96 such ladders in the complete detector. The four layers of the 
SVX are located at radii of 2.9 cm, 4.2 cm, 5.7 cm, and 7.9 cm from the beamline. 
Axial micro-strips, with a 60 jxm pitch on the three innermost layers and a 55 /Ltm 
pitch on the outermost layer, provide precision spatial measurements in the r tp- 
plane transverse to the beam. The geometric acceptance of the SVX is about 60% 
of the pp interactions, as it covers only ± 25 cm from the nominal interaction point, 
whereas the luminous region of the Tevatron beam has an RMS of ~30 cm along z. 

Because of radiation damage to the SVX readout chip, the performance of the 
SVX deteriorated over the course of the Run la data taking period. At the beginning 
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Fig. 8. (a) Isometric view of one of the SVX cylindrical modules, (b) Sketch of a ladder detector. 

of Run lb, the SVX was replaced with a new silicon vertex detector (SVX')El which 
is equipped with a radiation hard readout chip able to tolerate more than 1 MRad 
of radiation. The typical charge gain of the readout chip is around 21 mV/ f C, with 
typical values for the noise of around 1300 electrons (10.8 ADC counts), compared to 
about 2200 electrons for the original SVX detector. One improvement in noise level 
resulted from the fact that SVX' is AC coupled, in comparison to SVX which was 
DC coupled. This allowed the SVX' detector to be operated in double sample and 
hold mode, with one charge integration, resulting in a reduction in noise by a factor 
of a/2 compared to SVX which was operated in quadruple sample & hold mode, 
with two successive charge integrations for the determination of the outgoing signal. 

The SVX' position resolution is found to be better than 10 /im after alignment 
of the detector, using track data as can be seen in Fig. ||a). Here, the residual 
distribution of the distance of track intersection from reconstructed cluster centroids 
on a SVX' layer is plotted. The track impact parameter resolution is measured to be 
about (13 + 40/pt) ^mwithp t given in GeV/c. In Figure [|b) , the SVX track impact 
parameter resolution is shown versus 1/pt- Note the contribution from multiple 
scattering at lowp t , while the impact parameter resolution at highp t is dominated 
by the intrinsic detector resolution. The pt resolution of the CTC combined with the 
SVX is <r{p t )/pt = [(0.0066) 2 + (0.0009 pt) 2 ] 1 ' 2 , with p t measured in GeV/c, while, 
for CTC tracks alone, the resolution is a{p t )/p t = [(0.0066) 2 + (0.002 p t ) 2 } 1/2 - 

4. Collection of B Physics Data 

In this section, we describe the collection of datasets used for B physics analyses. We 
begin with the description of the B physics trigger data sets, followed by a summary 
of the selection requirements used to identify leptons, hadrons, and jets. At the end 
of this section, we describe the Monte Carlo simulation of CDF's B physics data. 

4.1. B physics triggers 

The total inelastic pp cross section at the Tevatron is about three orders of magni- 
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Fig. 9. (a) SVX' residual track position resolution, (b) Track impact parameter resolution 
versus 1/pt- 



tude larger than the b production cross section, putting certain requirements on the 
trigger system in terms of finding B decay products. In addition, the cross section 
for b quark production is a steeply falling cross section, as seen in Figure ^. It drops 
by almost two orders of magnitude between a b quark p t of about 8 GeV/c and 
25 GeV/c. In order to find B decay products in hadronic collisions, it is desirable 
to go as low as possible in the decay products transverse momentum, exploiting 
as much as possible of the steeply falling b cross section. Of course, the limiting 
factor is the bandwidth of CDF's data acquisition system. Throughout Run I, the 
CDF collaboration was able to maintain low p t -thresholds, without increasing the 
deadtime of the experiment during data taking. 

In Run I, all B physics triggers at CDF are based on leptons. Inclusive single 
lcpton triggers (e and fi) and dilepton triggers (dimuon and e\i) both exist. CDF 
uses a three-level trigger system to identify events to be written to tape. The first 
two levels are hardware based triggers. The Level 1 trigger uses information from 
detector subsystems, such as the muon chambers or the calorimeter, reducing the 
beam crossing rate from about 300 kHz to about 2 kHz. At Level 2, information 
from different detector subsystems is combined, such as a track stub in the CTC 
with hits in the muon chambers. This reduces the trigger rate to about 30 Hz. 
Level 3 is a software trigger based on the offline reconstruction code, optimized for 
computational speed, resulting in an event rate to tape of up to 10 Hz. 

4.1.1. Inclusive single lepton trigger 

The inclusive single lepton trigger identifies events with at least one electron or one 
muon, providing datasets enriched in events from semileptonic B decays. Inclusive 
electrons are selected at Level 1 by the presence of a single calorimeter tower above 
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Table 1 . Numbers of events recorded by the different B physics trigger streams for the Run la and 
Run lb data taking periods. 



Trigger Stream 


Run la 


Run lb 


Single electrons 




1.9 


10 6 


5.7 


10 6 


Single muons 




4.1 


10 5 


2.1 


10 6 


Dimuons 












< 2.8 GeV/c 2 




5.2 


10 5 


8.0 


10 5 


2.8 GeV/c 2 < m MM 


< 3.4 GeV/c 2 


2.2 


10 5 


1.0 


10 6 


3.4 GeV/c 2 < m MM 


< 4.1 GeV/c 2 


2.5 


10 5 


7.4 


10 5 


8.5 GeV/c 2 < m MM 


< 11.3 GeV/c 2 


8.1 


10 4 


2.8 


10 5 


m w > 4.0 GeV/c 2 




5.1 


10 5 


1.9 


10 6 


Electron-muons 




2.4 


10 5 


5.0 


10 5 



a threshold of 6-8 GeV depending on run conditions, while inclusive muons require 
the presence of a track in the CMU as well as the CMP. At Level 2, both of these 
triggers demand a charged track with p t > 7.5jGeV/c reconstructed in the r t/p-plane 
of the CTC by the central fast tracker (CFT)cZI, a hardware track processor, which 
uses fast timing information from the CTC as input. The momentum resolution 
of the CFT is cr(p t )/Pt = 3.5% with a high efficiency. In the case of the electron 
trigger, this track has to be matched to a cluster in the electromagnetic calorimeter 
with transverse energy E t > 8.0 GeV. In the case of the muon trigger, this track 
must be matched to a reconstructed track-segment in both the CMU and CMP. 
The trigger efficiency for a single lepton turns on at pt ~ 6 GeV/c, rises to about 
50% at a transverse momentum of ~8 GeV/c, and typically plateaus at over 90% 
at pt ~ 10 GeV /c depending on the exact trigger conditions. 

At Level 3, a computer farm is used to fully reconstruct the data, including 
three-dimensional track reconstruction in the CTC. However, the fast algorithm 
used for tracking is only efficient for particles with p t > 1.4 GeV/c. In the third 
level of the trigger, more stringent electron and muon selection criteria, similar to 
those described in Section 4.2, are applied. During Run I, about 7.5 ■ 10 6 electron 
trigger events and about 2.5 • 10 6 inclusive muon trigger events were recorded by 
CDF. Table |l| details the numbers of events recorded by the different trigger streams 
for the Run la and Run lb data taking periods. The numbers of events do not 
necessarily scale with the integrated luminosities of Run la and Run lb. This is 
due to different trigger thresholds and prescale factors, which accept only every 
second or third event triggered at high instantaneous luminosities, to prevent an 
increased deadtime. 

As an example, Figure |io|a) shows the transverse energy spectrum of electrons 
from the single electron triggers recorded in Run lb. The trigger turn-on is visible 
at -Et ~ 8 GeV followed by an exponentially falling energy spectrum. The shoulder 
around 40 GeV originates from a different physics process. These are electrons from 
W boson decays W — > ev, demonstrating how large the b cross section is compared 
to other physics processes like W production. The distribution in Fig. |ic|a) is 
obtained by identifying electrons with E t > 5 GeV in single electron trigger events. 
The enhancement of electrons at 5 GeV below the trigger turn-on results from 
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Fig. 10. Distribution of (a) transverse energy spectrum of inclusive single electrons and pt spectrum 
of (b) inclusive single muons and (c) muons from dimuon trigger data. The upper curve in (c) 
represents Run lb events while the lower distribution is obtained from Run la data. 



additional volunteer electrons in these events, which are not the trigger electrons. 
For comparison, Figure p~0|b) shows the pt-spectrum of inclusive single muons. 

4.1.2. Dimuon trigger 

The collection of dimuon trigger data requires two muon candidates be observed in 
the muon system at Level 1. The trigger efficiency for a muon at Level 1 rises from 
about 50% at p t = 1.6 GeV/c to 90% at Pt ~ 3.1 GeV/c with a plateau of -95%. 
As an illustration of the trigger turn-on, Figure pT|a) shows the Level 1 efficiency 
for CMU muons plotted versus 1/pt using J/tp — > /i/i and Z° — > /x/i data. The 
range-out at low momenta is predicted to occur for p t < 1.4 GeV/c. The second 
level trigger requires that at least one of the muon tracks is matched in ip to a 
track found by the CFT. The efficiency for finding a track with the CFT rises from 
50% at p t of about 1.9 (2.6) GeV/c to 90% at p t ~ 2.2 (3.1) GeV/c and reaches a 
plateau of ~94% (93%). The numbers given in brackets refer to the Run la settings 
where the dimuon trigger was operated with a slightly higher momentum threshold. 
Figure |ll|b) shows the Level 2 efficiency for CMU muons plotted versus 1/pt again 
using J /ip — > and Z° — > \x\x data. The Level 2 efficiency is displayed for positive 
and negative muons separately, indicating no charge dependence of the L2 trigger 
turn on. 

At Level 3, the events are again fully reconstructed, requiring two CTC tracks 
to be matched with two tracks in the muon chambers. For offline reconstruction 
the dimuon data are split into different samples, according to the dimuon invariant 
mass, as detailed in Table [l] which again shows the numbers of events separately for 
the Run la and Run lb data taking periods. The data with 2.8 < m w < 3.4 GeV/c 2 
are the source of CDF's J/ip — > events, while the data with 3.4 < m w < 

4.1 GeV/c 2 contain i(>(2S) — > events. Dimuon decays of the T resonances 
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Fig. 11. (a) Level 1 and (b) Level 2 trigger turn-on efficiency for CMU muons plotted versus 1/pt 
using J/ip —* fifi and Z° — > fifi data. The Level 2 efficiency is displayed for positive (top) and 
negative (bottom) muons separately. 



are contained in the sample with 8.5 < m MM < 11.3 GeV/c 2 . Events with to mm < 
2.8 GeV/c 2 are a source of double-semileptonic decays b — > c[iX followed by c — > 
/iJT, while the high-mass dimuon stream with m w > 4.0 GeV/c 2 provides events 
where the two muons originate from the decay of a b and b quark. 

As a further illustration, Figure p0|c) shows the pt-spectrum of muons from the 
dimuon trigger data. The upper curve in Fig. [To|c) represents Run lb events, while 
the lower distribution is obtained from Run la data. The higher muon momentum 
threshold in Run la is clearly visible and the trigger turn-on is softer compared to 
the single muon events (see Fig. [To|b). 

4.1.3. Electron-muon trigger 

CDF also triggered on events containing an electron and a muon. The trigger is a 
combination of the single electron and muon triggers. This sample has a principal 
electron E t threshold of about 5 GeV and a principal muon p t threshold of about 
3.0 GeV/c for Run la data and ~ 2.5 GeV/c for Run lb. Because of the higher 
lcpton threshold, the e/x-data comprise a significantly smaller dataset compared to 
the dimuon data as shown in Table [0. 



4.2. Lepton identification 



4.2.1. Electron identification 

The identification of electron candidates reconstructed after data collection uses in- 
formation from both the calorimeters and the tracking chambers. Here, we describe 
typical selection criteria used for electron identification at CDF. The longitudinal 
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shower profile must be consistent with an electron shower with a leakage energy from 
the CEM into the CHA of less than 4%, if one track is pointing to the calorimeter 
tower, or less than 10% if more than one track is pointing to the calorimeter tower. 
The lateral shower profile of the CEM cluster has to be consistent with that from 
test beam electrons. Additionally, a x 2 comparison of the CES shower profile with 
that of test beam electrons has to result in \ 2 < 10. For the association of a single 
high p t track with the calorimeter shower, based on the position matching at the 
CES plane, it is required that r|A<^| < 1.5 cm and |Azsin6>| < 3 cm. In addition, 
we demand the £/r of the electron candidate, reconstructed offline, to be greater 
than 6 GeV (5 GcV) for single electron (e/z) data. Usually electron candidates from 
photon conversion, are removed from the analysis by looking for oppositely charged 
tracks which have a small opening angle with the electron candidate. 

4.2.2. Muon identification 

The typical reconstruction of muon candidates at CDF is very similar for the single 
muon and dimuon datasets. Usually, a x 2 characterizing the separation between the 
track segment in the muon chamber and the extrapolated CTC track is computed, 
where the uncertainty in this \ 2 variable is dominated by the multiple scattering 
in the detector material. Typical selection requirements are \ 2 < 9 in the re- 
view (CMU, CMP, and CMX) and X 2 < 12 in the r-z view (CMU and CMX). A 
minimum energy deposition of 0.5 GeV in the hadron calorimeter is required for a 
muon candidate, as expected for a minimum ionizing track. The transverse muon 
momentum reconstructed offline is usually required to be greater than 7.5 GeV/c for 
inclusive muons. For the dimuon sample, this cut is typically placed at p t > 2 GeV/c 
for each muon candidate. 

4.3. Charged particle selection 

Charged particles are identified as tracks in CDF's central tracking chamber. To 
ensure a good track reconstruction, quality cuts requiring a minimum number of 
hits in the CTC are usually imposed. Typical requirements are at least five hits 
in two or more CTC axial superlayers and at least two hits in more than one 
CTC stereo superlayer. In addition, tracks are often required to be reconstructed 
in the SVX with hits in at least three out of the four silicon layers. In many 
analyses, a x 2 is calculated and defined as the increase in the track fit % 2 , when 
the SVX hits are included in the CTC track fit. This x 2 divided by the number 
of SVX hits usually has to be less than six to reject badly measured tracks. The 
track quality requirements are generally applied to hadrons as well as leptons. For 
optimal vertex resolution, hadron and lepton candidate tracks are often required to 
be reconstructed in the SVX detector. 

4.3.1. Hadron selection 

Usually, all reconstructed tracks are assigned the desired particle hypotheses, due 
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to the lack of a sophisticated particle identification at CDF. Applying cuts on the 
tracks transverse momenta are often the only means to reduce combinatorial back- 
ground. The minimum p t cut on tracks is usually set at 0.4 GeV/c. However, the 
specific ionization information dE/dx from the CTC provides a tt/K separation of 
about one a for tracks with p t greater than about 1.2 GcV/c. Therefore, dE/dx 
information from the CTC is sometimes used to help identify hadrons. Because of 
the large Landau tail of the ionization distribution, the 80% truncated mean of the 
measured charges from the CTC sense wires is taken as the best estimator of the 
track dE/dx. The probabilities, P(i), for a track to be consistent with the i — e, /i, 
ir, K, or p hypotheses arc then calculated using the measured dE/dx value and the 
predictions for the assumed particle hypotheses. A likelihood ratio, ih^ E ^ dx , for a 
track being, for example, a kaon is defined to be the ratio of P(K) divided by the 
sum of the probabilities of all particle hypotheses. Usually, this likelihood ratio is 
required to be greater than 0.01. 

The kinematic selection criteria used in a particular analysis are often optimized 
by maximizing the quantity Ns/\/Ns + Nb, where N$ is the predicted number of 
signal events based on Monte Carlo calculations (see Sec. 4.5) and Nb is the observed 
number of background events estimated e.g. from the signals sideband regions. 

4.4. Jet reconstruction 

Some analyses require the reconstruction of jets which arc usually formed from 
charged particles in B physics analyses, rather than from the more commonly used 
calorimeter clusters. A jet in a B event is often of low momentum and its energy 
is therefore not very precisely measured in the calorimeter. In addition, there is 
a difference between jets associated with electrons compared to jets containing a 
muon, since electrons deposit much more energy in the calorimeter than muons. If 
calorimeter clusters were used to find jets, electrons would be associated with jets 
much more often than muons for any minimum jet energy requirement. Therefore, 
jets are formed in B physics analyses from tracks using a cone clustering algorithm. 
Usually, tracks with p t > 1.0 GcV/c are used as jet seeds. If two seeds are within a 
cone with radius AR = \J (Arj) 2 + (Aip) 2 < 0.7, the momenta of the seeds are added 
together to form a new seed. After all possible seed merging, lower momentum 
tracks (usually 0.4 < p t < 1.0 GcV/c) that are within AR < 0.7 of a seed are 
added to form the final jet. A jet can, in principle, consist of a single track with 
Pt > 1 GeV/c. 

4.5. Simulation of heavy flavour production and decay 

Monte Carlo (MC) generation programs are used to simulate heavy flavour pro- 
duction in pp interactions. These programs are supplemented by phenomenological 
models of hadronization and decays of unstable particles. They are usually used 
in conjunction with programs that simulate the detector response to the final state 
particles. MC programs are necessary tools to calculate the geometrical and kine- 
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matic acceptance for observing heavy flavour events and are also used to estimate 
backgrounds for particular analyses. While considerable effort has gone into devel- 
oping and tuning simulations for e + e~ collisions, the state of the art is somewhat 
less well developed for the more complex high energy hadron-hadron collisions. 

At CDF, two types of Monte Carlo simulations are used. Calculations depending 
only on the production and decay of B hadrons employ a Monte Carlo generator 
that simulates only a single B hadron and its decay products. Situations which 
depend upon the fragmentation particles resulting from the hadronization of the 
b quark, as well as the "underlying event" particles from the pp scattering, use 
a full event generator like the PYTHIA program package!!. In the following, we 
describe both Monte Carlo generators as well as the CDF detector simulation. 

4.5.1. Simulation of a single B hadron 

The simulation of a single B hadron begins with a model of b quark production based 
on a next-to-leading order QCD calculationEj. This calculation employs the MRSDO 
parton distribution functior£3 with mj = 4.75 GeV, A = 215 MeV, and a renor- 
malization scale of /i = (j,q = \/p\ + m\ to model the kinematics of the initial state 
partons. Usually, b quarks are generated in the rapidity interval \y\\ < 1.0 with a 
minimum pt for the b quark chosen in a way to avoid any biases in the B meson kine- 
matic distributions after the application of the kinematic cuts used in an analysis. 
The b quarks are then fragmented into B hadrons, with no additional hadronization 
products, according to a model using the Peterson fragmentation functioned with a 
Peterson parameter of th — 0.006 (see Sec. 2.3). The bottom and charm hadrons are 
decayed into the various hnal states using branching ratios and decay kinematics 
governed by the world average masses and lifetimes of the invohzEd particlestH. The 
decay of bottom hadrons is accomplished using the QQ programed developed by the 
CLEO Collaboration, extended to include Bg mesons and b baryons. 

4.5.2. Monte Carlo simulation of the whole event using PYTHIA 

The PYTHIA Monte Carlo generatorfl is used in instances where more than just 
a single decaying B meson is required. PYTHIA simulates a complete pp interac- 
tion: The bb pair, the hadronization products, and the remaining beam fragments 
from the pp scattering ("underlying event"). PYTHIA exploits an improved string 
fragmentation model tuned to experimental data, mostly from high energy e + e~ 
collisions. PYTHIA generation at CDF uses most of the typical default parameters 
including the CTEQ2LE2I parton distribution functions. The b quarks are again frag- 
mented using the Peterson fragmentation modeled. However, the actual B decay 
performed by PYTHIA is suppressed and instead, the QQ programed is invoked. 

The PYTHIA generator is controlled by a series of parameters whose default 
values are adjusted to achieve good agreement with primarily high energy e + e~ 
data. Discrepancies between the "default" PYTHIA and CDF pp data are appar- 
ent, especially when considering particle production that does not originate from 



Collection of B Physics Data 21 



Table 2. The PYTHIA Monte Carlo parameters modified from their defaults to agree with 
B -> DtX data at CDF. 



Parameter 


Default 


Tuned 


Description 


MSTP(82) 


1 


3 


model of multiple interactions 


PARP(85) 


0.33 


1.0 


fraction of color-connected gg multiple interactions 


PARP(86) 


0.66 


1.0 


total fraction of gg multiple interactions 


MSTP(33) 


No 


Yes 


multiply cross sections by PARP(31) 


PARP(31) 


1.00 


1.66 


increase cross sections by 66% 


PARJ(21) 


0.36 


0.6 


<T fra S( P t) 


MSTJ(ll) 


4 


3 


use Peterson fragmentation for b, c 


PARJ(55) 


0.006 


0.0063 


Peterson fragmentation parameter ej 



the b hadronization but from the "underlying event" . The fidelity of the "default" 
PYTHIA generator is studied by comparing generated semileptonic B — » D£X 
Monte Carlo data after detector simulation to real CDF data. This comparison 
studied track multiplicities in AR and Atp intervals around the B — ► D£ direction 
and found the data to have a higher multiplicity of underlying event tracks than 
PYTHIA predicts. A good description of the charged particle multiplicities and 
Pt distributions is obtained by adjusting several PYTHIA parameters. The prop- 
erties of multiple interactions and beam remnants are controlled primarily through 
the multiple interaction cross section [PARP(31)], the model for their generation 
[MSTP(82)], the ratio of gg and qq multiple interactions [PARP(85,86)], and the 
width of the Gaussian px spread of particles produced in the breakup of color strings 
[PARJ(21)]. Once these parameters are adjusted to obtain agreement with the data 
away from the b jets, it is assumed that the underlying event is well modeled. In 
addition, the Peterson parameter PARJ(55) is modified so that the generated multi- 
plicity of tracks inside a cone of |Ai?| < 1 around the b quark matches the observed 
one. Table || lists the default and tuned values of the relevant PYTHIA parameters. 
More details of the tuning procedure may be found in Ref.Eil. 

4.5.3. Detector simulation 

The output of the MC generation is often simply passed through the standard fast 
simulation of the CDF detector which is based on parametrizations and simple 
models of the detector response determined from data or test beam measurements. 
The detector response is often parametrized as a function of the particle kinematics. 
The inclusive lepton or dilepton trigger usually introduces a strong kinematic bias 
on an analysis. This bias must be well modeled in the MC simulation to obtain the 
proper relative reconstruction efficiencies. Often, an empirical approach is taken 
rather than simulating the trigger directly. The trigger is either modeled by a 
simple parametrization of the CFT trigger, depending on the lepton p t , or by an 
error function parametrization of the ratio of the observed lepton p t distribution in 
the data compared to that generated by the simulation. For example, the efficiency 
parametrization of the single electron trigger follows the functional form 



22 B Lifetimes, Mixing and CP Violation at CDF 



Table 3. Comparison of important features of different experiments studying B physics. 





e+e" -> T(4S) -» BB 


e+e" -> Z° -> 66 


pp — » bbX 


Accelerator 


CESR, DORIS 


LEP, SLC 


Tevatron 


Detector 


ARGUS, CLEO 


ALEPH, DELPHI 


CDF, D0 






L3, OPAL, SLD 




o-(66) 


~ 1 nb 


~ 6 nb 


~50 (ib 


cr(bb) : cr(had) 


0.26 


0.22 


~ 0.001 


B°. B+ 


yes 


yes 


yes 


B%, B+, A° b 


no 


yes 


yes 


Boost < /3j > 


0.06 


6 


~2-4 


bb production 


both B at rest 


66 back-to-back 


66 not back-to-back 


Multiple events 


no 


no 


yes 


Trigger 


inclusive 


inclusive 


leptons only 



where A = 0.927, B x = 6.18 GeV/c, Ci = 4.20, B 2 = 7.48 GeV/c, C 2 = 0.504, and 
freq is the normal frequency function. After the simulation of the CDF detector, 
the same selection criteria applied to the data are usually imposed on the Monte 
Carlo events. 

5. Features of B Physics at a Hadron Collider: A Brief Tour 

In this section, we highlight some of the features of B physics at a hadron collider. 
We also make an attempt to describe, in an illustrative way, how CDF studies 
B decays. We will emphasize some of the tools used to find B decay products in 
hadronic collisions, omitting technical details which may be found later. 

To set the stage, we first compare current producers of B hadrons. Table | 
summarizes some of the important features of B physics experiments and the ac- 
celerators at which they operate. There are three main approaches to produce 
B hadrons. First, e + e~ — > T(4S f ) — > BB at the CLEO experiment, located at 
the CESR storage ring at Cornell and the DORIS storage ring at DESY, where 
the ARGUS experiment operated until 1993. Second, e + e~ — > Z° — > bb at the 
four experiments ALEPH, DELPHI, L3, and OPAL, located at the LEP storage 
ring at CERN, as well as the SLD detector at the SLC Collider at SLAC. Finally, 
pp — > bbX at the Tevatron where the CDF and D0 detectors are operated. As 
discussed earlier, the main motivation for studying B physics at a hadron collider 
is the large b quark cross section (see Tab. ||). 

Figure [l2] shows examples of typical B events at the T(4S) recorded with the 
ARGUS detector (left) and the CLEO experiment (right). At the T(45) resonance, 
only B°B° or B + B~ pairs are produced nearly at rest, resulting in a spherical event 
shape with an average charged particle multiplicity of about ten tracks. At the LEP 
or SLC accelerators, bb quark pairs are produced from the decay of the Z° boson 
where both quarks share half of the energy of the Z° resonance of 91.2 GeV. This 



results in two b jets being back-to-back as seen in Figure 13, where examples from 
the OPAL (left) and SLD (right) experiments are displayed. The average boost of 
B hadrons at the Z° resonance is f3j ~ 6. 

Figure |lj represents a typical B event from the Tevatron, shown on the left 
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Fig. 13. A typical B event at the Z° pole recorded with the OPAL detector (left) in the rip-view 
and the SLD experiment (right) in the rz-view. 
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hand side in the r ip-view of the CDF central drift chamber. No well defined jet 
structure is visible; the average multiplicity is about 40 charged tracks including 
tracks from the "underlying event" particles. It might appear challenging to find 
the B decay products in this quite messy environment of a hadronic collision. One 
way to extract B decays in a pp collision is illustrated on the right hand side of 
Fig. |l4|. The relatively long lifetime of a B hadron and a particle boost of about 
/?7 ~ 2-4, depending on the trigger which recorded the event, results in a B decay 
vertex which is clearly separated from the primary pp interaction vertex. In the 
example shown on the right hand side of Fig. Q the decay vertex of a Bg — > J/tp(p 
decay with J/tp — > fjrfiT and <p — > K + K~ is clearly visible and about 3.8 mm 
separated from the primary interaction vertex. 

In the following, we give a few more examples on how B decays are studied 
at CDF and which essential tools are used for B physics in hadron collisions. In 
Section 4.1, we already pointed out the importance of the trigger to select events 
containing B hadrons. Another important feature for B physics at a hadron collider 
is the good tracking capability such as that at CDF. The central tracking chamber 
together with the SVX provide excellent track momentum resolution which trans- 
lates into an excellent invariant mass resolution, as illustrated in Figure |l5|. Here, 
the dimuon invariant mass is displayed for muons from the dimuon trigger stream. 
A prominent J/tp peak is visible with about (243, 000 ± 540) J/tp signal candidates 
from the Run I data on low background. For this distribution, where both muons are 
reconstructed in the SVX, the mass resolution of the J/tp peak is about 16 MeV/c 2 . 

In addition to excellent tracking, superb vertexing is the other essential fea- 
ture of successful B physics studies at a hadron collider. This is demonstrated 
in Figure |l5|b) where the two muons of the J/tp signal candidates (light shaded 
area in Fig. |l5|a) are vertexed using tracking information from the SVX. The two- 
dimensional distance between the primary pp interaction vertex and the recon- 



Features of B Physics at a Hadron Collider: A Brief Tour 25 




Fig. 15. (a) Invariant mass distribution of oppositely charged muon pairs from the CDF dimuon 
trigger, (b) Decay length distribution of the signal events indicated as shaded area in (a) with the 
result of a fit superimposed (see text). 



structed dimuon vertex is plotted. This distribution shows several features: A 
prominent peak at zero decay length results from prompt J/ip candidates which are 
produced at the primary interaction vertex and constitute about 80% of all J/ip 
candidates. The width of this peak reveals information about the vertexing resolu- 
tion which is on average 40-50 /im, for this sample. At positive decay lengths, J/ip 
mesons from B hadron decays are described by an exponential slope. At a distance 
of about 100 (im from the primary interaction vertex, mainly J/ip's from B decays 
remain. There is also a small exponential slope at negative decay lengths where the 
particle seems to decay, before the point where it is produced. These events result 
from the combinatorial background underneath the J/ip signal. This is indicated 
by events from the J/ip sidebands (dark shaded regions in Fig. |l^a) which describe 
well the distribution at negative decay lengths as seen by the dark shaded area in 
Figure |l5|b) . 

In order to fully reconstruct B mesons, for example through the decay B + — > 
J/ipK + , the J/ip candidates from the signal region of Fig. |l5|a) are paired with 
another track in the event, assumed to be a kaon. A J/ipK + invariant mass dis- 
tribution, such as the one in Figure |l6|a), is observed with a B + signal on a large 
background. This background can be drastically reduced if a displaced B vertex is 
required. This is demonstrated in Fig. |l6|b) with a cut on the B decay length of 
greater than 100 /im. A clear B signal can now be found on a small background. 

Another tool in separating B decays in hadronic collisions, can be obtained from 
the fact that B decays are often isolated. This means there is usually not much 
track activity in the vicinity of a B hadron decay as illustrated in Fig. [l7|a). This 
fact can be exploited by a track based isolation quantity /. It is typically defined as 
the pt of the reconstructed B candidate divided by the scalar sum of the transverse 
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momenta of all tracks within a cone in 77 (£>-space around the B candidate direction, 
including the tracks from the B decay. A typical cone radius is AR < 1 and a typical 
cut value requires more than 50% of all momentum within the cone to be carried 
by the B candidate. The effect of the isolation cut is demonstrated in Fig. |l7|b) on 
a <fi — > K + K~ signal from a B® — ► Dg n + v decay with Dg — > cf)fi~D. The upper 
distribution shows the K + K~ invariant mass before a cut on / > 0.5 is applied to 
obtain the lower distribution. The combinatorial background is drastically reduced 
with a high efficiency on the <p signal. Finally, the significantly harder momentum 
spectrum of particles from b hadrons, compared to light hadrons, can be used to 
place high p t -cuts of several GeV/c on B candidates. This also results in a strong 
reduction of combinatorial background from light quark jet production. 

To summarize, the essential features to extract B hadron decays in hadronic 
collisions at CDF are excellent tracking in the CTC and SVX, the superb ver- 
texing capabilities of CDF's silicon vertex detector exploiting the long lifetimes of 
B hadrons by requiring decay vertices to be displaced from the primary interaction 
vertex, the harder momentum spectrum of particles from b hadrons, and the fact 
that B mesons are often isolated. 

6. Measurements of B Hadron Lifetimes 

CDF has measured the lifetimes of all weakly decaying B mesons (B°, B + , B$, 
and £?+) as well as the lifetime of the A" baryon. In this section, we discuss the 
B lifetime measurements at CDF. After a short introduction in the theory of B life- 
times, we describe experimental techniques of B lifetime measurements including 
the determination of the primary pp interaction vertex, the measurement of the 
B hadron decay length, and a description of the fitting procedure. We then review 
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Fig. 17. (a) Illustration of the isolation of B mesons in hadronic collisions, (b) Invariant K + K~ 
mass distribution from — > DJ fi + u decays with DJ — > 4>pu~v. The upper distribution shows the 
K + K~ invariant mass before an isolation requirement is applied resulting in the shaded K + K~ 
mass distribution. 



the lifetime measurements of the individual B species at CDF emphasizing crucial 
aspects of a particular measurement rather than describing an analysis in detail. 

6.1. Introduction to B hadron lifetimes 

The lifetimes of B hadrons are fundamental properties of these particles and are 
often needed to determine other Standard Modelcl quantities like the CKM matrix 
element \V c b\ or to measure the time dependence of BB oscillations. In the simple 
spectator model of a B hadron decay, where the b quark decays into a c quark 
emitting a W boson which couples to a Iv or qq pair, as illustrated in Figure [l^, 
the q quark within the B hadron only acts as a spectator and the b quark decays 
as a free particle. In this case, the B hadron lifetime would be given in analogy to 
the muon lifetime as 

r 192 7r J 

where T is a phase space factor. Here, we have neglected b — » u transitions resulting 
in a term with |T4tb| 2 which is small. In the spectator model, all hadrons containing 
a b quark would have the same lifetime. However, this picture does not hold for the 
prediction of charm hadron lifetimes which are measuredu to be 

t(D~) - 2.5t(L>°) - 2.5t(L>s ) - 5.0t(A~). (9) 

Possible causes for these lifetime differences originate from non-spectator effects, 
playing an important role in the decay of charm hadrons. These non-spectator de- 
cays include contributions like the W exchange process for neutral B mesons, shown 
in Figure |l^a), or the B + annihilation diagram displayed in Fig. |l9|b). Other 
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Fig. 18. Sketch of the spectator model of a B hadron decay. 
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Fig. 19. Illustration of (a) W exchange process, (b) B~*~ annihilation diagram as well as final 
state interference through (c) external W emission and (d) internal W emission. 



non-spectator mechanisms are caused by so-called final state Pauli interference ef- 
fects. An example is shown in Fig. 19 for the decay B + — > D°tt + which can occur 



through a so-called external W emission (Fig. |l9|c) or a so-called internal W emis- 
sion (Fig. |l9|d). There is a destructive interference between both diagrams which 
causes the B + decay width to be smaller than the B° width. This is one of the 
reasons the B + lifetime is predicted to be larger than the B° lifetime. Measure- 
ments of the lifetimes of individual B hadron species can therefore probe B decay 
mechanisms beyond the simple spectator model decay picture. However, among 
bottom hadrons, the lifetime differences are expected to be smaller than in the 
charm system due to the heavier bottom quark mass. 

In the past few years, the heavy quark expansion technique has been applied 
extensively to the calculations of inclusive decay rates of heavy hadrons, for both 
spectator and non-spectator decays. It provides quantitative predictions for life- 
time differences among the heavy hadrons. It is generally believed that a lifetime 
difference of order 5-10% should exist between the B + and B° meson. ReferencdiJ 
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predicts that the B + meson lifetime should be longer than the B° meson lifetime by 
about 5%. However, Referenced questions some assumptions made in Ref.EI and 
states that the sign of the deviation from unity cannot be predicted reliably. There 
is agreement that the models expect a much smaller difference of less than about 
1% between the B° and Bg lifetimes. 

Since this subject seems to be controversial, it might be best solved by precisely 
measuring all the B hadron lifetimes. Several direct measurements of B + and 
B° meson lifetimes have been performed by e + e~ experiments (see Ref.i for an 
overview) and by CDF. The precision of current measurements now approaches the 
level where the predicted small differences can be seen and improvements in these 
measurements will provide a strong test of B hadron decay mechanisms. 

6.2. Experimental techniques 

In general, two experimental methods have been employed to measure B hadron 
lifetimes. The first is based on the signed track impact parameter, the distance 
of closest approach of the track trajectory extrapolated to the B hadron produc- 
tion point. The average impact parameter is proportional to the lifetime of the 
B hadron. The advantage of using the impact parameter is that it is fairly insensi- 
tive to the boost of the B hadron: A B hadron with a large Lorentz boost will travel 
further, but the decay products will come out at a smaller angle leaving the impact 
parameter unchanged. To extract the B hadron lifetime, a Monte Carlo model is 
used to reproduce the observed impact parameter distribution as a function of the 
B lifetime. Impact parameter measurements typically usejeptons from semileptonic 
B decays. In fact, the first measurement of a 6 lifetimeo used the signed impact 
parameter of leptons. 

The second method for measuring B hadron lifetimes is based on the decay 
length, which is the distance from the B hadron production point to the B hadron 
decay point. The decay length L is related to the proper decay time t in the 
B restframe by the Lorentz boost f3"f as L = f3"f ct. Unlike the impact parameter 
method, it is necessary to know the boost value. In the case of a fully reconstructed 
B hadron decay, its boost value is determined as /?7 = pb/ttib, where ps is the 
B hadron momentum and mj is the B hadron mass. If the B hadron is only 
partially reconstructed, the boost value must be inferred using a Monte Carlo model. 

All B hadron lifetime results at CDF are based on the decay length measure- 
ment using fully reconstructed B hadrons as well as partially reconstructed B de- 
cays, which usually comprise higher statistics samples. In the following, we give 
an overview of how B lifetimes are measured at CDF, using the more complicated 
example of a partially reconstructed B decay. We first describe the Tevatron beam 
profile and the determination of the B hadron production point, the primary event 
vertex. Then, the B hadron decay length reconstruction using CDF's silicon vertex 
detector is subject of Section 6.2.3, while Sec. 6.2.4. gives a description of the life- 
time fit procedure used to extract the B hadron lifetime from the measured decay 
length. 
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Fig. 20. (a) Distribution of primary pp interaction vertices along the proton direction (z) for 
a typical data acquisition run during the 1994-1995 running period, (b) The two-dimensional 
distribution of the beam spot for a typical data acquisition run is shown in the upper two diagrams. 
The x- and y-projections, respectively, are shown in the lower two plots. 



6.2.1. Tevatron pp beam profile 

With a proton and antiproton bunch length of about 0.4 m, the primary pp interac- 
tion vertices are distributed along the beam direction according to approximately a 
Gaussian function, with a width of ~30 cm as shown in Fig. pO|a). Near the inter- 
action region, the p and p beams follow straight lines but could have an offset and 
slope with respect to the z-axis of the tracking detectors. The profile of the beam for 



a typical data acquisition run is shown in Fig. 20 o) . The deviation of reconstructed 
primary vertices from the calculated average beam position is plotted in the trans- 
verse plane. To ensure that the spread of the beam, instead of the resolution of the 
vertex fit, is the dominant contribution to the width of the observed distribution, 
only vertices with high track multiplicities are used. The upper two plots show the 
two-dimensional distribution of the beam spot for a typical Run I data acquisition 
run. The lower two plots show the x- and y-projections, respectively, with a fit to a 
Gaussian distribution superimposed. This shows that the transverse profile of the 
Tevatron beam was roughly Gaussian and circular with a width a of ~ 25 /zm in 
both the x- and y-directions. 

6.2.2. Primary event vertex 

The B hadron lifetime measurements reported in this article are based on measuring 
the distance between the primary pp event vertex and the secondary B decay vertex 
in the transverse plane exploiting CDF's silicon vertex detector. First, the z-position 
of the primary interaction vertex is identified, using the tracks reconstructed in the 
VTX detector and then the transverse position of the primary event vertex is deter- 
mined, using the run average beam position. The tracks reconstructed in the VTX 



Measurements of B Hadron Lifetimes 31 



detector, when projected back to the beam axis, determine the longitudinal location 
of the primary interaction with an accuracy of about 0.2 cm along the beam direc- 
tion. On average during Run I, the number of reconstructed interaction vertices in 
a given event follows a Poisson distribution with a mean of about 2.5. For B hadron 
lifetime measurements, the z-location of the primary event vertex is usually deter- 
mined by choosing the pp interaction vertex, recorded by the VTX which is closest 
to the intercept of the B candidate with the bcamlinc. Usually, the z-coordinates 
of all tracks from the B decay are required to be within five cm of the z-position 
of this primary vertex, in order to exclude tracks from other pp interactions in the 
event. 

The transverse position of the primary event vertex is determined using the 
average beam position through the detector, together with the knowledge of the 
longitudinal primary vertex position from the VTX. The average beam position 
is calculated offline for each data acquisition run. It is found that the average 
beam trajectory is stable over the period of which a given pp beam is stored in 
the Tevatron Collider. For B lifetime measurements, usually only events from data 
runs with a sufficiently large number of collected events are considered in order to 
allow a good determination of the run averaged beamline. 

In B lifetime analyses, the primary vertex is usually not measured on an event- 
by-event basis because the presence of a second b quark decay in the event, coupled 
with the low multiplicity in e.g. semileptonic B decays can lead to a systematic 
bias in the lifetime determination. The algorithm that determines the primary 
event vertex on an event-by-event basis, calculates first the z-position of the pri- 
mary pp interaction vertex from the VTX in the same way described above. The 
transverse position of the primary vertex is then determined for each event by a 
weighted fit of all SVX tracks, with a z-coordinate within 5 cm of the z-vertex 
position obtained from the VTX. At first, all SVX tracks are fit to originate from 
a common vertex. Tracks which have large impact parameters, with respect to this 
vertex, are then removed and the fit is repeated in an iterative process until no more 
tracks are removed. At least five tracks must remain for a successful primary vertex 
reconstruction to occur. The uncertainty in the fitted primary vertex coordinates 
transverse to the beam direction ranges from about 10-35 /im depending on the 
number of tracks used in the fit and the event topology. 



6.2.3. Decay length measurement 

In the following, we describe the reconstruction of the B hadron decay length on the 
more complicated example of a partially reconstructed B decay. We shall point out 
the difference to the simpler case of a fully reconstructed B hadron. As our example, 
we choose a semileptonic decay B — ► D£v which would represent a lifetime mea- 
surement using, for example, B° — > D~ ' l + v or B% — ► Dgl + v decays. A schematic 



representation of such a semileptonic B — > Dlv decay is displayed in Figure 21. 

Usually, the D candidate is searched for in a cone in 77 93-space around the lepton 
candidate and fully reconstructed. The tracks forming the D candidate are then refit 
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Fig. 21. Schematic representation of a semileptonic B — > D£y decay. 

with a common vertex constraint referred to as the tertiary vertex and marked as Vd 
in Fig. The secondary vertex, where the B decays to a lepton and a D (referred 
to as Vb), is obtained by simultaneously intersecting the trajectory of the lepton 
track with the flight path of the D candidate. If the D meson is fully reconstructed, 
which is usually the case, the D flight path and thus the B decay vertex are known. 
The confidence level (C.L.) of the combined vertex fit is typically required to be 
greater than 1%. Furthermore, it is often required that the reconstructed D decay 
vertex Vd be positively displaced from the primary vertex, as projected along the 
direction of the Dt momentum. 

The decay length L^ y is then defined as the displacement of the secondary vertex 
Vb from the primary event vertex Vp, measured in the plane transverse to the beam 
axis and projected onto the transverse momentum of the D£ system (see Fig. |2l| ): 



t b _ (VB-Vp)-Pt(Dl) nm 



L^y is a signed variable which can be negative for the configuration, where the 
particle seems to decay before the point where it is produced. The B meson decay 
time is given by 

,„s r m(B) 
Pt{B) 



where m(B) is the mass of the B hadron. In the case of a fully reconstructed 
B hadron, Pt(B) is known and ct (B) is used as input to fit directly for the B life- 
time. In the case of a semileptonic B decay, where the B meson is not fully recon- 
structed, a "pseudo-proper decay length" is usually defined as 



m{B) 
Pt(D£)' 



A = £f y ^- (12) 



Using single lepton trigger data, the reconstructed B decay length has a typical 
uncertainty of ~ 50-60 /im, including the contribution from the finite size of the 
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Fig. 22. Examples of normalized i<"-factor distributions pt(D s i + )/pt(Bg), for Bg — > D s 
Monte Carlo decays with (a) DJ" — » </<7r~ and (b) Dg — > (f>fi~u. 



primary event vertex. In addition, the ratio if of the observed momentum to the 
true momentum is introduced as 

K= MbJ> (13) 

to correct between the reconstructed pt(D£) and the unknown pt(B) in the data. 
The B meson decay time is then given as 



The correction between p t (D£) and pt{B) is done statistically by smearing an ex- 
ponential decay distribution with a Monte Carlo distribution of the correction fac- 
tor K, when extracting ct(B) from the pseudo-proper decay length in the lifetime 
fit procedure as described in the next Section 6.2.4. 

The Zf-distribution is obtained from D£ combinations which originate from a 
Monte Carlo simulation of the semileptonic B decay of interest. As an example, 
the -fT-distribution is shown in Figure ^ for the Bg lifetime measurement using 
the semileptonic decay Bg — > Dg£ + v which is further described in Sec. 6.5. Fig- 



ures 22 1) and |22p) show the iC-distribution for the cases where the Dg meson is 
reconstructed in the Dg — * <fiir~ and Dg — > cf)fi~D decay modes, respectively. These 
isT-distributions have mean values of 0.86 and 0.77 with RMS values of 0.10 and 0.12 
for Dg — > cjm~ and Dg — > (j)pL~v, respectively. The i^-distribution is approximately 
constant as a function of p t (D£) in the range of interest, which is typically 15 GeV/c 
to 25 GeV/c for single lepton trigger data. 

To ensure a precise B lifetime determination, usually only B candidates are con- 
sidered for which the pseudo-proper decay length is measured with an uncertainty of 
less than 0.1 cm. In addition, the D candidates are often required to have a proper 
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decay length measured between Vb and Vd of less than 0.1 cm and an uncertainty 
on this proper decay length of less than 0.1 cm. These requirements reject poorly 
measured decays and reduce random track combinations. 

6.2.4. Description of lifetime fit procedure 

We now describe the fit procedure used to extract the B lifetime from the pseudo- 
proper decay length distribution. We continue to use our example of a semileptonic 
B — > Dtv decay. As input to the B hadron lifetime fit, usually, a signal sample is 
defined by a mass window around the respective D mass peak. To model the pseudo- 
proper decay length distribution of the background events contained in the signal 
sample, a background sample is defined consisting of events from the D sidebands. 
See Fig. |l5|a) for a better understanding of the definition of signal events (light 
shaded area) and sideband samples (dark shaded area), shown here for a J/tp — > 
invariant mass distribution. 
The pseudo-proper decay length distribution obtained from the signal sample is 
then fit using an unbinned maximum log-likelihood method. Both the B lifetime, 
denoted as cr below, and the background shape are determined in a simultaneous 
fit using the signal and background samples. Thus the likelihood function £ is a 
combination of two parts 

N s N B 

£ = n if**?** + o- - fs^H, } ■ n 

* 3 

where Ns and Nb are the number of events in the signal and background samples. 
Here, / s i g is the ratio of D signal events to the total number of events in the signal 
sample. To constrain / s j g , which is usually a free fit parameter, to the number of 
D signal events (/sig) with uncertainty er s ; g , obtained from the D mass distribu- 
tion, an additional x 2 term x 2 = (/ s ; g — (/ s i g )) 2 /er 2 ig is factored into the likelihood 
function defined above. 

The signal probability function F s i s consists of a normalized decay exponential 
function convoluted with a Gaussian resolution function Q and is smeared with a 
normalized if-distribution 'H(K) 



r sig (x) = J dKH(K) 



K Kx 

— exp{ } ® Q{\i\ x, sai 

CT CT 



(16) 



Here, A; is the pseudo-proper decay length measured for event i with uncertainty 
<7j and x is the true pseudo-proper decay length. The symbol "<g>" denotes a con- 
volution and Q(Xi \ x, sai) is the Gaussian distribution given by 

g(Xi\ x, sat) = 1 -= exp (~ ^z ^ V (17) 



Because of systematic uncertainties in the overall scale of the decay length uncer- 
tainties, which are estimated on an event-by-event basis, a scale factor, s, is often 
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introduced. It is usually a free parameter in the B lifetime fit. The integration over 
the momentum ratio K is approximated by a finite sum 

J dKH(K) -» AKH{Ki), (18) 

i 

where the sum is taken over bin i of a histogrammed distribution TL{Ki) with bin 
width AK as shown e.g. in Figure ^2| 

The background probability function is typically parametrized by a Gaus- 
sian centered at zero, a negative exponential tail, and a positive decay exponential 
to characterize the contribution of heavy flavour decays in the background sample: 

- f+ - f-)G(Xi\x,s(Xi) + -±-exp{-— } ® g(X i \x,sa i ) + 

X 

exp{- — } ® G(\i\ x, so-j). (19) 

Here, f± are the fractions of positive and negative lifetime backgrounds and A± 
are the effective lifetimes of those backgrounds. The parameters usually allowed to 
float in the fit are the B lifetime, / S i g , X±, f±, and the overall scale factor s. 

For a fully reconstructed B decay, the smearing with a normalized X-distribution 
TL{K) does not apply but the convolution with the Gaussian resolution function Q 
is essential to obtain the true pseudo-proper decay length x. In the case of a fully 
reconstructed B decay, the background is usually also described by a Gaussian 
centered at zero, a negative exponential tail, and a positive decay exponential. 

6.3. Measurement of average B hadron lifetime 

The measurement of the average B hadron lifetime is based on B — ► J/ipX decays 
where the J/ip is reconstructed through its decay intojL + /i~. This measurement 
constitutes the first publication of a B lifetime at CDFE3 using only about half of 
the Run la data. This analysis demonstrated in 1993 the capabilities of the_newly 
installed silicon vertex detector. The final average B hadron measurementcS uses 
a subset of the J/ip data sample already shown in Figure p"5|a). Only the Run lb 
data and only J /tp candidates, where both muons are matched to a CFT track at 
Level 2, are used. Since the precision of the measurement is limited by systematics, 
including more data would not improve the result. 

This analysis uses the two-dimensional decay length L xy of the J/ip vertex and 
corrects for the difference between the boost of the J/ip meson and the B hadron, us- 
ing a Monte Carlo model which compared the background subtracted p t -distribution 
of the data with the J/ip momentum of the MC simulation (see Sec. 6.2.3.). Care is 
taken in this analysis to avoid non-Gaussian components in the errors on the decay 
length distribution. To ensure that the J/ip vertex is well measured and to have a 
good understanding of the a of the vertex fit, strict track and vertex quality cuts 
are applied. These required, for example, both muons to have associated clusters in 
all four layers of the SVX or the calculated uncertainty on the decay length <JL xy to 
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be less than 150 /im. This leaves 67,800 pairs of oppositely charged muons in the 
J/ip signal region and 7,900 pairs in the combined sidebands, which are used to de- 
termine the decay length distribution of the background underneath the J/ijj signal. 

The pseudo-proper decay length distribution for the signal region, with the result 
of the fit superimposed, is already shown in Fig. p~5|b) . The dark shaded area 
shows the contribution from background, where the shape is determined from the 
sidebands and the magnitude is derived by normalizing the number of sideband 
events within Poisson fluctuations to the same range in invariant mass as used for 
the signal sample. The light shaded region shows the contribution due to adding the 
exponential distribution from B decay to the background. The remaining unshaded 
region, which constitutes more than 80% of the J/tp sample, shows the contribution 
from prompt J /if) mesons produced at the primary interaction vertex. 

The measured inclusive B lifetime, which is the average over all B hadrons 
produced in pp collisions at ^/s = 1.8 TeV weighted by their production cross 
section, branching ratios and detection efficiencies, is 

(t(B)) = (1.533 ± 0.015 ±°;°3i) PS- (20) 

The first uncertainty is statistical and the second is systematic. This is the standard 
way to present results where two uncertainties are quoted in this paper. 

6.4. Measurement of B + and B° lifetime 



6.4.1. B + and B lifetimes with fully reconstructed B mesons 

The analysis principle for the B lifetime measurement using fully reconstructed B + 
and B° mcsonsE3 is as follows. All Run I dimuons forming a J ft]) candidate, as shown 
in Figure |l5|a), are used. B + mesons are reconstructed in the decay modes J/iJjK + , 
J/?pK* + , i(j(2S)K + , and ip{2S)K* + , while B° mesons are reconstructed in the 
decay modes J/<ipK%, J/iPK*°, ^(2S)K%, and ^(25")^*° with ip(2S) -> /i+ \i~ . The 
kaons are reconstructed using the decay channels K*° — > K + ir~, K* + — » Kgir~, 
and Kg — > 7r + 7r~. K* candidates are accepted if their invariant Kir mass is within 
±80 MeV/c 2 of the world average K* massi. Since swapping the assignment of the 
kaon and pion masses to the two tracks forming a K*° candidate can also result in a 
K*° candidate within the ±80 MeV/c 2 mass window, the K*° combination closest 
to the world average mass value is chosen. 

The vertex and mass constraint J ftp and ip(2S) candidates are fit with the 
kaon candidates to originate from a common vertex, yielding the two-dimensional 
B meson decay length L^ y . Together with the measured B transverse momentum p t , 
L^ y is used to obtain the proper time distributions shown in Fig. ^3| for (a) charged 
and (b) neutral B candidates. The bottom plots represent the backgrounds as 
obtained by fitting the B sideband regions to a Gaussian with exponential tails. 
Using this background shape, an unbinned maximum log-likelihood fit of the signal, 
assumed to be an exponential convoluted with a Gaussian, is performed. The 
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Fig. 23. The proper decay length distribution for fully reconstructed (a) B + mesons and (b) 
B° mesons. The upper (lower) histogram shows the peak (sideband) region distribution. The su- 
perimposed curves are the contributions from the signal, the background and their sum determined 
by the log-likelihood fit. 



following lifetimes are obtained: 

t{B+) = (1.68 ± 0.07 ± 0.02) ps, 
t(B°) = (1.58 ± 0.09 ± 0.02) ps, 
t(B + )/t(B°) = 1.06 ± 0.07 ± 0.02. (21) 

When calculating the uncertainty on the lifetime ratio, the correlated systematic 
errors are properly taken into account. 

The exclusive B lifetime measurement using fully reconstructed B decays is 
statistics dominated. One way to increase the number of B candidates is not to 
fully reconstruct the B meson. This is done in the semi-exclusive analysis described 
in the next section. 

6.4.2. B + and B° lifetimes with partially reconstructed B mesons 

The B lifetime analysis using partially reconstructed B mesonsEll exploits the semi- 
leptonic decays B — > D^*'£i/X and follows the description given in Sec. 6.2.3 and 
6.2.4, as illustrated in Fig. [2l]. The analysis starts with the single lepton trigger 
data and searches for charm mesons in a cone around the trigger electron or muon. 
DM meson candidates are reconstructed through their decay modes: 

(a) D° — > K + ir~ , where the D° is not from a D*~ decay, 

(b) D*~ -> D°7T-, D° -> K+tt-, 

(c) D*~ D°n-, D° — ^ K+n-n+n-, and 

(d) D*~ — > D°7r~, D° — * K + ir~Tr°, where the 7r° is not reconstructed. 

The charm signals for the decay modes (a)-(d) above can be seen in Fig. |5|a)-d), 
respectively. Figure |2~i|a) shows the K + n~ mass spectrum while in Fig. ^Jb)-d) the 
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Fig. 24. Charm signals reconstructed in the B lifetime measurement using partially reconstructed 
B — > D(*}£vX decays with (a) D° — » iC+7r~, where the D° is not from a D*~ decay, (b) D*~ — » 
£) 7r-, D° -» ^T+tt-, (c) D*" -> D°7T-, D° -» E'+tt-tt+tt- , and (d) D*~ -> D°7r-, D° -* 
K + Tr~~n°. The shaded histograms show wrong-sign or D°£+ combinations. 



mass difference Am between the measured invariant masses of the D°ir~ and D° 
candidates is displayed. The mass peak is broadened for mode (d) because of the 
missing 7r° meson. Note, the charm signals in Fig. ^ are quite clean and rather 
competitive with signals found at e + e~ colliders. The number of events in the 
signal regions and the estimated background fractions are summarized in Table ^[ 

The D° and D*~ candidates are then intersected with the lepton to find the 
B decay vertex. Since the B meson is not fully reconstructed, a /?7 correction has 
to be applied to scale from the observed £)Wf momentum to pt(B) as described 
in Sec. 6.2.3. The obtained lifetime distributions from D°£ + and D*~£ + are used 



Table 4. Summary of the definition of signal samples, numbers of charm candidates and estimated 
background fractions. 



B Mode 


D° mode 


D° mass range 


Am range 


Events 


Background fraction 






[GeV/c 2 ] 


[GeV/c 2 ] 






l+D° 


K+TT- 


1.84 - 1.88 


Not D*~ 


5198 


0.526 ± 0.018 


e+D*~ 


K+n- 


1.83 - 1.90 


0.144 - 0.147 


935 


0.086 ±0.011 


e+D*~ 


K ^ 7T~ TY^ 7T~ 


1.84 - 1.88 


0.144 - 0.147 


1166 


0.183 ± 0.015 


t+D*- 


K+n-n° 


1.50 - 1.70 


< 0.155 


2858 


0.366 ± 0.016 
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Fig. 25. Pseudo-proper decay length distribution of (a) D°£ + and (b) D*~l + candidates where 
the three D° decay modes are combined. The curves show the result of the combined fit: The 
B° component (dot-dashed curves), the B + component (dashed curves), and the background 
component (dotted curves). 



to determine the individual B + and B lifetimes. A D°£ + combination usually 
originates from a charged B meson while a D*~£ + candidate comes from a B°. This 
simple picture is complicated by the existence of D** states which are the source of 
D° (£>*-) mesons originating from a decay B° -> D**-£+X, £)**- -> D Q X (B+ 
jj**0£+X, D**° — > D*~X). This cross talk from D** resonances, decomposed using 
Monte Carlo, is an important aspect of this analysistll. A combined lifetime fit to 
the pseudo-proper decay length distributions, using an unbinned maximum log- 
likelihood method, as described in Sec. 6.2.4, yields the results shown in Fig. ^5|. 
The contamination of the wrong B species is only at the 10-15% level. It is indicated 
in Fig. ^5], in addition to the respective signal and background contributions. 

The final lifetimes of the B + and B° mesons, using their partially reconstructed 
semileptonic decays B + — > D°£ + vX and B° — » D*~£ + vX, are 

t(B+) = (1.637 ± 0.058 i°;^)ps, 
r(B°) = (1.474 ± 0.039 iO; 05?)PS. 
t(B+)/t(B°) = 1.110 ±0.056 IS, (22) 

where the main systematic error is from the estimate of the B meson momentum. 
6.4.3. B° and B + lifetime summary 

We combine the B lifetime measurement using fully reconstructed decays, described 
in Sec. 6.4.1 and given in Eq. ([jl|), with the CDF measurement using partially 



reconstructed B mesons, described in Sec. 6.4.2 and given in Eq. (22), to derive the 
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Fig. 26. Comparison of the CDF B° 
ratio measurements (right) with other 
as of July 1998. 



lifetime measurements (left) and the T{B Jr )/r(B ) lifetime 
experiments, compiled by the LEP B lifetime working group 



following CDF average B lifetimes: 

T (B+) = (1.661 ± 0.052) ps, 
t(B°) = (1.513 ±0.053) ps, 
t(B+)/t(B°) = 1.091 ±0.050, (23) 

where the uncertainties include both statistical and systematic effects. There ex- 
ists a small (about 0.01 ps) correlation in the systematic effects between the two 
measurements. This is taken into account in combining the results. The ratio of 
the two B meson lifetimes differs from unity by about 9%, or almostTwo standard 
deviations. This agrees with the small difference predicted by theoryE3. 

As an example, a comparison of the CDF B° lifetime measurements with other 
experiments, compiled by the LEP B lifetime working group in July 1998 for the 
29th International Conference on High Energy Physics, Vancouver, Canada, can 
be found in Figure |2^ on the left hand side. A comparison of the CDF B lifetime 
ratio measurements compared to other experiments is presented in Figure ^ on the 
right hand side. These comparisons show that the CDF B lifetime measurements are 
very competitive with results from the Z° pole at LEP and SLC. The current world 
average B° and B + lifetimes as of July 1998, determined by the LEP B lifetime 
working group, taking correlated systematic uncertainties into account are: 

t(B+) = (1.67 ±0.03) ps, 
t(B°) = (1.57 ±0.03) ps, 
t(B+)/t(B°) = 1.07 ±0.03. (24) 
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plus a polynomial background, (b) Bg proper decay length distribution for events selected within 
±0.05 GeV/c 2 of the fitted mass. The results of the lifetime fit are shown: Signal (dashed 
line), background (dotted line), and the sum of the two (solid line). 



6.5. Measurement of Bg lifetime 

CDF has measured the Bg lifetime using fully and partially reconstructed Bg decays. 
A search for a lifetime difference AT /T in the Bg system is described in Section 6.5.3. 

6.5.1. Bg lifetime with fully reconstructed Bg mesons 

In this exclusive Bg lifetime measurement!!, the Bg candidates are reconstructed 
in the decay chain Bg — > J/tpcp, with J/tp — > and <p — > K + K~ . The dimuon 

data sample used in the lifetime measurement with fully reconstructed B and 
B + mesons (see Sec. 6.4.1) is also the starting point to reconstruct J/tp mesons in 
this analysis. Once a J/tp is found, <j> — > K + K~ candidates are searched by selecting 
oppositely charged track pairs assigning each track the kaon mass. If the mass of 
the <f> candidate lies within ±10 MeV/c 2 of the world average (p massta, all four 
tracks are constrained to come from a common vertex. The J/tp(p invariant mass 
spectrum is shown in Fig. |27|a). A fit to a Gaussian plus a second order polynomial 
background yields a signal of (58 ± 12) events. 

Due to the more limited number of Bg candidates in this analysis, a maximum 
log-likelihood fit is simultaneously performed to the invariant mass and the proper 
decay length distributions, in order to better constrain the number of Bg signal 
events. A plot of the proper decay length spectrum is given in Fig. |27]b) for candi- 
dates within ±0.05 GeV/c 2 of the fitted Bg mass. The Bg lifetime is measured in 
the exclusive mode J/tpcp to be 

t(B°s) = (1.34 t° f 9 ±0.05) ps. (25) 
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Table 5. Summary of the four Dg decay modes: The estimated number N(Dg) of Dg signal 
events, the number JV ev t of events in the signal samples, and the fitted Bg lifetimes cr(Bg), where 
the errors shown are statistical only. 



Dg Decay Mode N(D~) N cvt ct(B°) 



(pTV 


220 ± 21 


350 


418 


+43 

-39 ^ m 


K*°K~ 


125 ± 20 


820 


411 


+73 

-66 A™ 


K%K- 


33 ±8 


146 


397 


till Mm 


(f>li~& 


205 ± 38 


635 


399 





Since this Bg lifetime measurement, using fully reconstructed Dg decays, is statis- 
tics dominated, using partially reconstructed Dg decays improves the statistical 
precision as described in the following section. 

6.5.2. Dg lifetime with partially reconstructed Bg mesons 

The Bg lifetime measurement with partially reconstructed Bg mesonsil uses the 
semileptonic decay Bg — > Dg£ + v and is very similar to the B + and D° lifetime 
analysis using partially reconstructed B decays (see Sec. 6.4.2). The Dg candidates 
are reconstructed in the decay modes 

(a) ->0tt-, <f>^K+K-, 

(b) D s -> K*°K-, K* Q -> K+TT-, 

(C) Dg -> K S K~, K° s -» 7T+7T-, 

(d) Dg -> fa~L>, <j> -> K+K~, 

For the first three decay modes the analysis starts with a single lepton trigger 
data set, while the semileptonic Dg decay mode is based on a dimuon data sample 
obtained with a trigger requirement of m(/x/i) < 2.8 GeV/c 2 . Dg candidates are 
searched for in a cone around the lepton and then intersected with the lepton to find 
the Bg decay vertex (see also Fig. |2l]). Since the Bg meson is not fully reconstructed, 
its c t cannot be directly determined and a f3"f correction is applied to scale from the 
Dg£ + momentum to pt(Bg). Figure ^ shows the Dg invariant mass distributions 
for (a) 07T-, (b) K*°K-, (c) K° S K- , and (d) K+K~ from Dg -> fa- v. The dots 
with error bars are for right-sign Dg£ + combinations while the shaded histograms 
show the corresponding wrong-sign Dg i" distributions which show no signals. The 
numbers of Dg signal events are compiled in Table |^. 

One of the crucial aspects of this analysis is that the reconstructions of the Dg 
decay modes into K*°K~ and KgK~ suffer from reflections of D~ — * K*°tt~ and 
D~ — > Kgir~ , respectively, where the ir~ is incorrectly assigned the kaon mass. We 
will discuss this reflection from D~ and the determination of the true number of 
events from the Dg decay with the example of the Dg — » K*°K~ mode. The effect 
of this K-tt misassignment can be seen in Figure p9|. Events from a B — > D~ IvX 
Monte Carlo simulation with D~ — > K*°ir~ yield an invariant mass distribution 
indicated by the shape of the shaded area in Fig. ^2$p ) if they are reconstructed as 
Dg — ► DglvX with Dg — > K*°K~, misinterpreting the tt~ as K~. A significant 
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portion of this D reflection lies at the D s mass peak. Two methods are used to 
determine the D~ reflection from the data. 

The first method performs a simultaneous fit to the K*°K~ and K*°ir~ invariant 
mass distributions, where the K*°ir~ mass distribution is created by switching the 
mass assignment on the K~ track to a pion. Figure ^9|a) shows the K*°K~ mass, 
while the corresponding K*°ir~ mass is displayed in Fig. 29 3). Each distribution 
is described by a Gaussian for the corresponding D~ and Dg signal, as shown 
in Figures ^c) and p9|d), plus a linear lineshape to parametrize the combinatorial 
background. The shape of the corresponding D~ or D$ reflection, as obtained from 
a Monte Carlo simulation, is also included in the fit as displayed in Fig. p9|c) and 
d) as the shaded areas. The two mass distributions are fit simultaneously with the 
number of events in the Gaussian Dg {D~ ) signal set equal to the number of events 
in the corresponding D~ (Dg) reflection. The fit result is shown in Figure p9|a)-d). 

The second method for determining the amount of D~ reflection in the Dg signal 
sample exploits the difference between the D~ lifetime [t(D~) — (1.057±0.0015) ps] 
and the D s lifetime [t(D s ) = (0.467 ±0.0017) ps]. Fitting for the Dg lifetime, the 
exponential describing the Dg signal is replaced by the sum of two exponentials, one 
with the -Dg lifetime and one with the D~ lifetime. The Dg and D~ lifetimes are 
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Fig. 29. (a) Mass distribution for candidates in the Dg —> K*°K~ decay mode, (b) Mass 
distribution if these candidates are assumed to be D~ — > K*°tt~ . (c) Distribution of the Dg — » 
K*°K~ signal and the reflection from D~ — » K*°it~ (shaded area) as obtained from Monte 
Carlo simulations. Normalizations are determined from the simultaneous fit described in the 
text, (d) Mass distribution of the corresponding D~ — » K*°n~ signal and the reflection from 
Dg —> K*°K~ (shaded area), (e) pseudo-proper decay length distribution for Dg —* 4>ft~ 
with the results of the fit superimposed. The dashed line is the Bg signal contribution, while the 
shaded area represents the contribution from the combinatorial background, (f) Pseudo-proper 
decay length distribution for the background sample with the fit result superimposed. 



fixed to their nominal valuesa and the relative fractions of Dg and D~ are allowed 
to float in the fit. The weighted average of Dg events from both methods yields a 
Dg signal fraction of (62 ± 10)% for the Dg — > K*°K~ decay. Both methods are 
also used to calculate the number of Dg events and the contribution from the D~ 
reflection in the Dg — » KgK~ decay mode with the result displayed in Table 0. 

The Bg lifetime is determined for each of the four Dg decay channels individ- 
ually, with the fit results and their statistical uncertainties shown in Table ^[ As 
an example, the pseudo-proper decay length distribution of the Dg — > cj)ir~ signal 
sample with the result of the fit superimposed is shown in Figure |29|e) . The dashed 
line represents the Dg signal contribution, while the shaded area shows the sum of 
the background probability function over the events in the signal sample. The same 
distribution of the background sample is displayed in Figure p^f). The combined 
Dg lifetime from all four Dg decay modes is determined from a simultaneous fit 
to be 

r(Dg°) = (1.36 ±0.09±° : °05)PS- (26) 

This result is currently the world's best measurement of the Dg lifetime from a 
single experiment. 

Using the CDF average D° lifetime r(D°) = (1.513 ± 0.053) ps (see Sec. 6.4.3), 
we determine the Dg/D° lifetime ratio to be 0.899 ± 0.072, taking correlated sys- 
tematic uncertainties into account. However, ignoring the correlated systematic 



Measurements of B Hadron Lifetimes 45 



uncertainties increases the error on the Bg/B° lifetime ratio only to ±0.077, since 
t(Bs)/t(B°) is dominated by the statistical error on the Bg lifetime measurement. 

6.5.3. Determination of AT /T 

In the Standard Model, the Bg meson exists in two CP-conjugate states, \Bg) = \bs) 
and \Bg) = \bs). The two mass eigenstates of the Bg meson, Bg and Bg (H = 
'heavy' and L = 'light'), are not CP eigenstates but are mixtures of the two CP- 
conjugate quark states: 

\Bf) = 1/V2(\B° S ) - |P°) ) and |P£) = 1/V5( \B° S ) + \B° S ) )■ (27) 

The mass and lifetime differences between the Bg and Bg can be defined as 

Am = m H -m L , Ar = T L -r H , and r = 1/2 (T H + T L ), (28) 

where mn,L and Tjj.l denote the mass and decay width of Bg and Bg. Please 
note, we have defined both Am and Ar to be positive quantities. In this case, the 
heavy state is the long-lived state or the CP odd state, while the light state is the 
short-lived state or the CP even state in analogy to the neutral kaon system. Un- 
like the B° meson, thcwidth difference in the Bg system is expected to be largeE3. 
Theoretical estimatesoO predict AT/T to be on the order of 10-20%. In the Bg sys- 
tem the ratio Am/ AT is related to the ratio of the Cabibbo-Kobayashi-Maskawafli3 
matrix elements |V c t,V cs |/|Vt s Vt(,|, which is qu ite w ell known, and depends only on 
QCD corrections within the Standard ModeEilc3. A measurement of Ar would 
therefore imply a determination of Am and thus a way to infer the existence of Bg 
meson oscillations, which will ultimately determine the ratio of the CKM matrix 
elements |Vtd|/|Vt s |. The current limit on Am in the Bg system!, Arris > 10-2 ps -1 
at 95% C.L., gives rise to expect a significant lifetime difference in the Bg system. 

An illustrative wayJxi better understand the relation between Am and Ar in 
neutral meson systemsc2l is shown in Figure |30[ Let's at first concentrate on the 
well known B°B° system displayed in Fig. |30|a). Here, the B° state is illustrated as 
its mass excitation curve, a Breit-Wigner lineshape, with the ir-axis representing an 
energy in ps -1 . The width of the Breit-Wigner curve is given by the reciprocal of 
the B° lifetime 1/T(B°) ~ 1.6 ps. The second curve represents the heavy B° state 
which is separated from the light B° state by a mass difference Am ~ 0.47 ps -1 
which is well measured from B°B° oscillations. The lifetime difference Ar is very 
small in the B° system resulting in both curves appearing with the same width. 

The classic system to illustrate Am and Ar is the K°K° system, displayed in 
Fig. |30|b). Here, a short lived state exists, the Kg, shown as the broad curve with 
a width of l/T(K s ) ~ 0.09 ns, and a long lived state exists, the K^, shown as the 
narrow curve with a width of 1/T(K%) - 52 ns. The lifetime of the K° L is about 
500 times longer than the Kg lifetime. The mass difference Am between the states 
is about 5.3 ns -1 . Please note, the unit here is in ns -1 compared to ps -1 for the 
B°B° system. This means particle antiparticle oscillations are much slower in the 
K° system compared to the B° system. 



46 B Lifetimes, Mixing and CP Violation at CDF 




Fig. 30. Illustration of the relation between Am and Ar in (a) the B°B° system, (b) the K°K° 
system, (c) the B^B® system, and (d) the D°D° system. 



The system of interest, in our case, is the BgBg system with a large mass 
difference between the Bg and Bg state. We choose e.g. Am = 15 ps -1 as drawn 
in Fig. |30|c). In analogy to the kaon system, the Bg is the shorter lived state, while 
the Bg is the longer lived state. In Fig. |30|c) a lifetime difference of AT/T ~ 0.3 is 
assumed resulting in the Bg state to appear slightly broader than the Bg . Finally, 
in the D°D° system, shown in Fig. |30|d), the mass difference Am is very small in the 
Standard Model as is Ar. This results in completely overlapping mass excitation 
curves for both states which appear as one curve with a width of 0.4 ps. 

It is assumed that Bg mesons are produced as an equal mixturecil of Bg and 
Bg . One way to search for Ar in the Bg system is to describe the Bg meson decay 
length distribution by a function of the form 



T{t) = l/2(r ff e~ r " t + r L e- rtt ) with T L , H = T ± AT/2, (29) 
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rather than by just one exponential lifetime re _rt . Clearly, the integral over time 
/ dtT(t) returns the same number of Bg and Bg mesons as produced. This would 
be the right functional form in a search for a lifetime difference AT/T, if a com- 
pletely inclusive sample of Bg mesons is used which is allowed to decay into any 
final state. However, for a semileptonic final state, the semileptonic decay widths 
are the same for Bg and Bg , since the reason for a lifetime difference originates 
from hadronic Bg decays. This can be easily seen in the following way: In an ex- 
cellent approximation, the semileptonic Bg decay is flavour specific, which means 
that a £ + can only originate from a \bs) state. Since the fraction of \bs) states is the 
same for Bg and Bg, the semileptonic widths are the same for Bg and Bg, thus 
the semileptonic branching ratios are different for Bg and Bg. The correct func- 
tional form to describe the Bg meson decay length distribution from semileptonic 
Bg decays, assuming a lifetime difference AT/T, is therefore: 

AT 1 AT 

^) = r H r L /r.(e- r « t + e- r - t ) with r L , ff = r±— = r-(i±-— ). (30) 

The parameter AT/T is the parameter fit for. 

In the case of a lifetime difference Ar 7^ 0, the total decay width T = 1/2 • 
(Th + Tl) and the mean Bg lifetime r m (_Bg) obtained from a fit assuming a single 
Bg lifetime, are no longer reciprocal to each other but follow the relation 

1 i + (AL)2 

'-<*>- F-ir{$ji- < 31 > 

The relation in Eq. ( |3l| ) is incorporated into the likelihood fitting function for AT/T, 
and the mean Bg lifetime is fixed to the world average B° lifetime, since both 
lifetimes are expected to agree!!!! within 1%. 

The fit returns AT/T — 0.34 where the given error is statistical only. 

This indicates that with the current statistics CDF is not sensitive to a Bg lifetime 
difference. Based on this fit result, the normalized likelihood is integrated as a 
function of AT/T and the 95% confidence level limit is found at 

AT/T < 0.83 (95% C.L.). (32) 

Using a theoretical value of AT/ Am = (5.6 ± 2.6) • 10~ 3 from Ref.0 and setting 
T m (Bg) to the world average B Q lifctimea, an upper limit on the Bg mixing frequency 
of Arris < 96 ps _1 (95% C.L.) can be determined within the Standard Model. 
Including the dependence on AT / Am and T m (Bg) into the limit, we obtain 

6.6. Measurement of A° lifetime 

The analysis principle for the A° lifetime measuremen is very similar to the B°, 
B + , and Bg lifetime analyses using partially reconstructed decays. The A° baryon 
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Fig. 31. (a) Invariant mass distribution of pK~ ir + for right-sign (points with error bars) and 
wrong-sign events (shaded area), (b) Pseudo-proper decay length distribution of events from the 
AjJ signal region. The points with error bars are data, the solid line is the fit result, the shaded 
area is the signal distribution, and the dashed line is the background contribution. 



is reconstructed through the semileptonic decay A® — * A+£~DX, with the subse- 
quent decay A+ — > pK~ir + . The analysis again uses the single lcpton trigger data 
searching for A+ candidates in a cone around the lepton. The A+ candidates are 
intersected with the lepton to find the A° decay vertex. A signal of (197 ± 25) A+ 
events is obtained as shown in Fig. |3l|a), where the pK~n + invariant mass distri- 
bution for right-sign A+£~ combinations is plotted. The shaded histogram shows 
the wrong-sign A+£ + distribution. 

The pseudo-proper decay length distribution of events from the A^ signal region 
is shown in Fig. |3l|b). The points with error bars are data, the solid line is the 
fit result, the shaded area is the signal distribution, and the dashed line is the 
background contribution. Using these events, the A° lifetime is determined to be 

r(Ag) = (1.32 ± 0.15 ± 0.07) ps. (34) 

The CDF A° lifetime result is competitive with the LEP measurements in 
precisionu. Using the CDF average B° lifetime t(B°) = (1.513 ± 0.053) ps (see 
Sec. 6.4.3), we determine theJifetime ratio r(Ag)/r(£f°) = 0.87 ± 0.11. Theory 
favours the value for this ratic£3 to be in the range 0.9-1.0 in good agreement with 
the CDF measurement. 

6.7. Measurement of £?+ lifetime 

Discussing the measurement of the Bf lifetime involves summarizing the recent 
discovery of the B+ meson at CDFo'EJ. The £>+ meson is the lowest-mass bound 
state of a bottom antiquark and a charm quark: \B+) = \bc). This pseudoscalar 
ground state has non-zero flavour and no strong or electromagnetic decays. It is 
the only state with two different heavy quarks where each can decay weakly. As 
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Fig. 32. (a) Sketch of quark periodic table, (b) Histogram of the J/tp£ mass distribution that 
compares the signal and background determined in the likelihood fits to the combined data for 
J/ipe and J/ip/i. Note that the mass bins, indicated by tick marks at the top, vary in width. The 
inset shows the log-likelihood function versus the number of HJ mesons. 



seen in Figure 32 1), the \bc) state is the last meson state predicted by the Standard 
Modela, which had not been discovered until early 1998. Figure ^a) shows a quark 
periodic table with all possible combinations of meson states which can be obtained 
from the five quarks u, d, s, c, and b. 

Non-relativistic potential models predict a masscll around 6.2-6.3 GeV/c 2 . 
In these models, the c and b quarks are tightly bound in a very compact system and 
have a rich spectroscopy of excited states. The predicted lifetimes of the B+ meson 
is in the range 0.4-1.4 ps. Because of the wide range of predictions, a lifetime 
measurement is a test of the different assumptions made in the various calculations. 

We attempt here a simple educated guess of the mass and its lifetime. We 
would estimate the B+ mass as 

m(B+) ~ m(b) + m(c) - 4.75 GeV/c 2 + 1.55 GeV/c 2 = 6.3 GeV/c 2 , (35) 

where we have approximated the b and c quark masses with half of the mass of 
the T(15) and the J/ip, respectively. To estimate the B+ lifetime, we expect three 
major contributions to the B+ decay width: T(B+) ~ Tb + T c + Tb c - Decays of the 
b quark b — > cW + with the c quark as a spectator, leading to final states such as 
J/tpTT or J/ip£v; c quark decays c — ► sW + , with the b as spectator, leading to final 
states such as Bgir or Bglv; and cb — > W + annihilation, leading to final states like 
D K , tv t or multiple pions. In the simplest view, the c and b decay like free quarks 
with no annihilation contribution, and the B+ lifetime would be: 

t(b+) = i/r(s+) - [r 6 + rj- 1 - [1/1.65 P s + 1/0.47 ps]- 1 - 0.4 ps , (36) 

where is approximated with the B + lifetime, while we use r(Dj) to estimate T c . 
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Table 6. Summary of signal and background events. 

3.35 < m(J/tpe) < 11.0 GeV/c 2 
J/ipe Events J/tp fj, Events 

False electrons 
Undetected conversions 
False muons 
BB background 
Total background (predicted) 

(from fit) 
Predicted N(B+ -> Jipeu)/N(B^ 
e and fi signal (derived from fit) 

Total signal (fitted parameter) 

Signal + Background 21.2 ±4.3 19.0 ±3.5 

Candidates 23 14 

Probability for null hypothesis 0.63 X 10~ 6 



4.2 ±0.4 
2.1 ± 1.7 

11.4 ±2.4 

2.3 ±0.9 1.44 ±0.25 
8.6 ±2.0 12.8 ±2.4 
9.2 ±2.0 10.6 ±2.3 

Ji/>tv) 0.58 ± 0.04 

12.0+^ 8A+H 



6.7.1. Discovery of B^ meson at CDF 

The meson is reconstructed via its semileptonic decay B+ — > J/ip£ + i/X (£ = 
e, fi). In this analysis, advantage is taken of a clean J/tp — > [i + signal (see Fig. |l5| 
and Sec. 6.4.1) and excellent lepton identification at CDF. To reduce backgrounds 
from prompt J/ ip production, tri-lepton vertices displaced from the primary vertex 
are searched for by requiring the J/ip£ pseudo-proper decay length to be greater 
than 60 /im. 23 B+ — » J /ipev candidates and 14 B+ — * J /ipfif candidates are 
found. Because of the missing neutrino the B^ cannot be fully reconstructed and 
no B+ mass peak can be obtained. The signal is therefore found as an excess of 
events over expected backgrounds. Gre at care is taken to correctly determine the 
different backgrounds in this analysis^. 

Significant backgrounds in the B+ candidates come from false leptons. Hadrons 
reaching the muon detectors without being absorbed, hadrons that decay in flight 
into a muon in advance of entering the muon detectors, and hadrons falsely identified 
as electrons. Background from photon conversions 7 — > e + e~ arises when one 
member of the pair remains undetected and the other accidentally intersects the 
J/tp decay point. Background from BB decays arises when a J/ip from a B decay 
and a lepton from a semileptonic decay of the R appear to originate from the same 
decay point. A number of other backgrounds^ are found to be negligible. Table ^| 
summarizes the results of the background calculation and of a simultaneous fit for 
the muon and electron channels to the mass spectrum over the region between 
3.35 and 11 GeV/c 2 . Figure |32]b) shows the mass spectra for the combined J/tp£ 
candidates, the combined backgrounds and the fitted contribution from — > 
J jiplv decays. The fitted number of U+ events is 20.4^5 5. 

To test the significance of this result, a number of Monte Carlo trials is gener- 
ated with the statistical properties of the backgrounds, but with no contribution 
from B+ mesons. These are subjected to the same fitting procedure to determine 
contributions consistent with the signal distribution arising from background fluc- 
tuations. The probability of obtaining a yield of 20.4 events or more is 0.63 x 10~ 6 , 
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Fig. 33. (a) Pseudo-proper decay length distribution for the combined J/ipe and J/i/'A' data 
with the fit result and the contributions from signal and background overlaid. The inset shows 
the log-likelihood function versus ct. (b) Measured value of the <r X BR ratio shown as point 
with one standard deviation contour and plotted at the measured value for r(B^~). The shaded 
regions represent theoretical predictions and their uncertainty bands for two different values of the 
B^t semileptonic decay width r s j. 



equivalent to a 4.8 standard deviation effect. 

To check the stability of the B+ signal, the assumed mass is varied from 
5.52 to 7.52 GeV/c 2 . The signal template for each value of m(B+) and the back- 
ground mass distributions are used to fit the mass spectrum for the data. The 
magnitude of the B+ signal is stable over the range of theoretical predictions for 
m(B+). The minimum in the log-likelihood function versus mass yields a B+ mass 
of m{B+) = (6.40 ± 0.39 ± 0.13) GeV/c 2 . 

The lifetime of the B+ meson is obtained from the decay length distribution 
using only events with 4.0 < m(J/ip£) < 6.0 GeV /c 2 and relaxing the requirement 
on the pseudo-proper decay length from > 60 /mi to greater than —100 /^m. This 
yields a sample of 71 events, 42 J/ipe and 29 J/V'M candidates. An unbinncd 
maximum log-likelihood fit yields 

t(B+) = (0.46 ±g:i| ±0.03)ps. (37) 

The data, together with the signal and background distributions, are shown in 
Fig. ||a). 

From the 20.4 £?+ events and a sample of 290 B + — > J/%l)K + events selected 
with the same requirements as the B^ candidates, the B+ production cross section 
times the B+ — > J/ip£ + i/ branching fraction a x BR(B+ — > J/tpl + v) is obtained, 
relative to that for the topologically similar decay B + — > J/tpK + . Many systematic 
uncertainties cancel in the ratio, while Monte Carlo calculations yielded the values 
for the efficiencies that do not cancel. The detection efficiency for — > J/i/ji+v 
depends on ct because of the cut on the pseudo-proper decay length at > 60 /im. 
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Fig. 34. Summary of CDF B hadron lifetime results. 



We therefore quote a separate systematic uncertainty from the lifetime uncertainty. 
We find 



a(B c ) x BR(B C -> J/ip£v) 
a{B) x BR(B -> J/tpK) 



0.132±°;°& (stat.) ± 0.031 (syst.) +»%? (cr) (38) 



for and B + with transverse momenta pt > 6.0 GeV/c and rapidities \y\ < 
1.0. Figure |33|b) compares phenomenological predictions with the measurements of 
t(I?+) and the a x BR ratio. Within experimental and theoretical uncertainties, 
they are consistent. 



6.8. B Lifetimes: Summary 

A summary of the B lifetime measurements at CDF is given in Figure |34|. As 
we have seen, CDF's B lifetime results are very competitive with the LEP and 
SLC measurements, where a precision of a few percent is reached. Although CDF's 
measurement of the B + / B° lifetime ratio appears to be different from unity by 
almost two standard deviations, the precision is still not yet sufficient to distinguish 
between theoretical approaches. The A" lifetime lies closer to the B° lifetime at 
CDF (t(A^]/t(B°) = 0.87 ± 0.11) compared to the measurements of the LEP 
experimental Finally, the B+ lifetime, measured for the first time at CDF, is 
clearly much shorter than the other B hadron lifetimes and closer to the D° or 
Dg lifetime. 
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7. Measurements of B°B° Oscillations 

Oscillations between particles and antjnarticles were predicted^ in 1955 and ob- 
served in the system of neutral kaonsEll in 1956. Oscillations in the B°B° meson 
system were observed for the first time by the ARGUS collaboration!!!! in 1987. This 
discovery signaled a large top quark mass and gave rise to new prospects for the 
observation of CP violation in B meson decays (see Sec. 8). In this section, we 
describe measurements of the time dependence of B°B° oscillations at CDF. After 
a short introduction in B flavour oscillations, we review different approaches for 
measuring the time dependence of B°B° oscillations at CDF. 

7.1. Introduction to B°B° oscillations 

As already outlined in Sec. 6.5.3, the system of neutral B mesons, B° and B°, can 
be described in terms of states with well defined mass and lifetime, Bh and Bl 
(H = 'heavy' and L = 'light'), 

\B°) = 1/V2 ( \Bl) + \B° H ) ) \B°) = 1/V2 ( \B Q L ) \B%) ). (39) 

Here, we consider the effects of CP violation to be small compared to the expected 
mixing effects. The difference in mass Am between both mass eigenstates leads to 
a time dependent phase difference between their wave functions. The probability 
that an initially pure B° state can be observed as a B at proper time t is given by 

"Pmix = Prob(B° ->B°,t) = 1/2 re"" (1 - cos Am t), (40) 

while the probability that it decays as B° is given by: 

7> un mix = Prob{B° -> B°,t) = 1/2 Te" rt (1 + cos Amt). (41) 

Figure [35] shows, on the left side, the time evolution of BB oscillations displaying 
the unmixed (solid) and mixed (dashed) contributions for two different oscillation 
frequencies Am. The sum of V m ix and "P un mix is just the exponential particle decay 
re _r * and is displayed as a dotted line in Fig. 

If the B meson decay time t cannot be measured, the time integrated probability 
which is usually referred to as x can be obtained as 

V ■ x 2 
x = = - (42) 

Pmix ^unmix 2(1 4" X ) 

Here, the mixing parameter x — Am/T is introduced. It describes the oscillation 
period relative to the B meson lifetime tb = h/T. 

If both neutral B mesons, B a and Bg, are produced, the time integrated and 
flavour averaged mixing parameter x is defined as 

X = fdXd + fsXs, (43) 

where fd and /g are the fractions of b hadrons that are produced as B° and 
Bg mesons, respectively. 



54 B Lifetimes, Mixing and CP Violation at CDF 




Fig. 35. Left: Time evolution of BB oscillations displaying the unmixed (solid) and mixed 
(dashed) contribution as well as the sum of both (dotted) for two different oscillation frequencies. 
Right: Box diagram describing B°B° oscillations. 



In the Standard Modelta, B°B° mixing occurs via second order weak processes, 
as displayed in Fig. |3^ on the right hand side. The mass difference Amd can be 
determined within the Standard Model by computing the electroweak box diagram, 
where the dominant contribution is through top quark exchange: 

Am d = ^| m B (f 2 B B B )rj QCD m\ F(^-) \V t * b V td \ 2 . (44) 

D7T m W 

Here, Gf is the Fermi coupling constant, raj the B meson mass, f B the weak 
B decay constant, Bb the bag parameter of the B meson, ?7qcd are QCD corrections 
which are in the order of one, m t is the top quark mass, Vtt and Vtd are the two 
CKM matrix elements involved, and F(z) is a slowly varying function of the top 
quark mass and the W boson mass mw- 



1 9_ 

4 + 4(1 -z) 2(l-2) 2 2(l-z) 3 



= 7 + 77T— X-oTT— ^-oTT— W (45) 



In the case of the B$ meson the respective values of B B o , f B o , m B c i , and »7qcd = Vb° 
have to be used, and Vtd is replaced by Vt s - 

A measurement of Am^ or Amg would in principle determine the Cabibbo- 
Kobayashi-Maskawa matrix elements Vtd or V ts but theoretical uncertainties con- 
nected with the poor knowledge of the B meson weak decay constant fs and the 
bag parameter Bb limit the direct extraction of the CKM matrix elements from 
measurements of Am^ and Arris- However, in the ratio Am^/Ams several of the 
theoretical uncertainties cancel 

Aw d _ m B o T] B a f B gB B o \V td \ 2 

2 R_„ IT/. 12- ^ 



Am s m B o r] B o fl B B o \V t 
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and Amj/Ams will ultimately determine one of the legs of the unitarity triangle 
as further discussed in Sec. 8.1. 

7.1.1. Ingredients of B°B° mixing measurements 

In general, a measurement of the time dependence of B°B° oscillations requires the 
knowledge of the proper decay time t of the B meson, the flavour of the B meson 
at production, and the flavour of the B meson at decay in order to determine 
whether the B° has oscillated. The last two items require identifying the flavour 
of a B meson. They are the subject of the next section 7.1.2, where we discuss 
B flavour tagging. As outlined in Sec. 6.2.3, the B meson decay time can be 
obtained from a measurement of the distance L^ y between the primary interaction 
vertex, where the B meson is produced, and its secondary decay vertex. The decay 
time t is related to the decay distance L^ y by 



ct(B) = ^=2£^. 

h y pt(B) 



If the B meson is not fully reconstructed, a /?7 correction is applied to scale from 
the only partially measured B momentum to the unknown B momentum. 

From Eq. ( fl7| ) the uncertainty on the decay time can be calculated (in units of 
the B lifetime tb): 



Hi 

tb 



/AL^\ 2 / t Ap 



where L° —Pt/mB • ct~b- (48) 



The proper time resolution a t depends on the uncertainty to infer the decay length 
from the primary to the B decay vertex and on the B momentum resolution which 
is dominated by the magnitude of the (3"f correction. Note, the latter uncertainty 
scales with i/ra, while the vertexing resolution only adds an uncertainty constant 
with ct. 



7.1.2. B flavour tagging 

B flavour tagging refers to the task of determining the flavour of a B meson either at 
production or at decay. Several methods of B flavour tagging exist. Some methods 
identify the flavour of the other B hadron produced in the initial collision along 
with the B meson of interest. Since the dominant b quark production mechanism 
produces bb pairs, the flavours of both B hadrons are assumed to be opposite at time 
of production. There are three common methods of opposite side flavour tagging. 
One method, called "lepton tagging" , looks for a lepton from the semileptonic decay 
of the other B hadron in the event. The charge of this lepton is correlated with the 
flavour of the B hadron: An l~ comes from a b — > cl~vX transition, while an £ + 
originates from a b quark. Second, the charge of a from the subsequent charm 
decay c — > sX is also correlated to the B flavour: A K~ results from the decay chain 
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b — > c — > s while a K + signal a b flavour. Searching for a charged kaon from the other 
B hadron decay is usually referred to as "kaon tagging" . Third, a method called "jet- 
charge tagging" exploits the fact that the sign of the momentum weighted sum of 
the particle charges of the opposite side b jet is the same as the charge of the b quark 
producing this jet. Finally, the flavour of a B meson can also be tagged on the same 
side as the B meson of interest, by exploiting correlations of the B flavour with the 
charge of particles produced in association with the B meson. Such correlations are 
expected to arise from B quark hadronization and from B** decays. We call this 
method "same side tagging" . CDF has results on opposite side lepton tagging, jet- 
charge tagging, and same side tagging but has not studied kaon tagging due to the 
lack of particle identification at CDF. These tagging methods are described in more 
detail in the following sections where they are applied in different B a B° mixing 
measurements. 

The figure of merit used to compare the "tagging power" of different flavour tags 
is the so-called effective tagging efficiency el? 2 , where e is the efficiency of how often 
a flavour tag is applicable. The dilution T> is defined as the number of correctly 
tagged events Nr minus the number of incorrectly identified events Nw divided by 
the sum 

N R + N W 1 ' 

To quantify the tagging power with the expression "dilution" can be misleading. A 
flavour tag which always tags correctly has a dilution of one, while a flavour tag 
giving the correct tag 50% of the time has a dilution of zero. This means, a tagging 
algorithm with a large dilution is desirable, while a small dilution characterizes a 
less powerful tagging method. The dilution is also related to the probability pn 
that the flavour tag is correct and to the mistag probability pw = 1 — Pr that the 
flavour tag is incorrect 

V = 2p R -l = l-2p w . (50) 

To illustrate the statistical significance of the product eT> 2 , we discuss an asym- 
metry measurement on a data sample with N events, where the flavour tagging 
method identifies whether the event is of type a or type b. Type a and type b 
could for example be "mixed" and "unmixed" decays of a neutral B meson. The 
measured asymmetry *4 mea s is 

Ancas ~ N a + N b ' (51) 
where N a and Nb are the number of events that are tagged as type a and type b, 
respectively. The true asymmetry *4 t ruc is 

where TV* and N b are the true number of events of type a and type b in the sample. 
The efficiency is simply 

Ng+Nb 

£ = WTM- (53) 
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It follows, from the definition of the dilution, that the true asymmetry and the 
measured asymmetry are related by Atmc = -4 mcas /2?, and the statistical error on 
the true asymmetry is 

h - V 2 A 2 

where N is the total number of events in the sample N — iV* + N£ . The product 
eV 2 N is the effective statistics of the data sample, that is, it is the equivalent 
number of perfectly tagged events. The statistical power of different flavour tagging 
methods varies as eT> 2 . 

7.1.3. Outline of B°B° mixing measurements at CDF 

All B°B° mixing measurements at CDF are based on lepton data samples where 
the charge of the lepton determines the B flavour at decay assuming it originates 
from a semilcptonic B decay. The B flavour at production is determined by an 
opposite side lepton tag, a jet-charge tag, or a same side tag in the various analyses 
described in the next sections. The proper time at decay is determined from the 
B decay vertex, inferred from partially reconstructed B mesons or inclusive vertices. 
To illustrate the effects discussed above, we consider a mixing measurement where 
an opposite side lepton identifies the B flavour at production. 

We start with a pure sample of B° mesons and assume that the lepton tag is 
always correct. In this case, an event with an opposite-sign lepton pair signals an 
unmixed event, while a like-sign lepton pair indicates a mixed event. The probabil- 
ities for an opposite-sign event Vos an d a like-sign event Vls are directly related 
to the mixing probabilities 

V LS {t)=V w i*{t) and Vos(t) = V unm Ut) (55) 
which are defined as 

Puamix/mixW = 1/2 r e ~ r ' (1 ± cos Am t). (56) 
To illustrate this behavior, we plot in Fig. |36|a) the asymmetry 

"J^0 CJ 1 — L S 

Anix = ^ ^ — = cos Am t (57) 

Vos + V LS 

which in this case is a pure cosine like oscillation. We have chosen Am = 5. Next, 
we introduce a vertexing resolution which smears the decay time measurement and 



effectively reduces the amplitude of the oscillation as shown in Fig. 363). The effect 
of introducing a momentum resolution, for example through a fC-factor distribution 
7~t(K) in case of a partially reconstructed decay, is displayed in Fig. |36|c). Since the 
uncertainty on the /3~f resolution scales with proper time t, as seen in Eq. (E^), the 
oscillation damps with proper time. In the case of a Gaussian vertexing resolution Q 
and a momentum resolution expressed through TL{K), the mixing probability would 
be modified as 

V mix (t)=f dKH(K) 1/2 KTe- Krt ' (1- cos Am Kt') ®g(t'\t,a). (58) 
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Fig. 36. Illustration of various effects on the asymmetry *4 m i x : (a) Perfect resolution, (b) 
vertex resolution, (c) momentum resolution, (d) mistag probability, (e) vertex and momentum 
resolution plus mistag probability, and (f) vertex and momentum resolution plus mistag including 
combinatorial background. 



T'mix would be modified accordingly. 

If the lepton tag in our example does not always tag correctly but with a mistag 
probability pwt we would measure the following like-sign and opposite-sign fractions 

VLs(t) = (1-pw) Vmix(t) + PW 7\nmi x (i) 
'Posit) = PW ^mix(t) + (1 PW ) ' unmix 

(i). (59) 

The effect of a mistag probability py/ is shown in Fig. |3^d), while the combined 
effects of vertexing plus momentum resolution, together with a mistag probability 
are displayed in Fig. ^6)s). In a real measurement there will be background such as 
combinatorial background under a charm signal. We define V^s an d Vqs as Pls 
and Vos as in Eq. (|5^) and obtain 

V LS (t) = (l-f bs )Vi s s (t) + h s hsV hB (t) 

Vos(t) = (l-.f bg )V s ^ s (t) + f hg (l-f LS )V hg (t). (60) 

Here, fh s is the fraction of background in a given sample, fog is the fraction of 
like-sign events in the background, while Vbg is the probability function to describe 
the behavior of the background versus proper time t. The effect of a background 
on the asymmetry, in addition to a vertexing and momentum resolution, as well as 
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v Candidate track P m 

Fig. 37. (a) A simplified picture of b quark fragmentation, (b) Schematic drawing of the 
momentum vectors determining the p T t el of a same side tagging track candidate. 



a mistag rate are shown in Fig. |36|f) and exhibit a distribution as can be expected 
in a real measurement. Some measurements demonstrate the time dependence of 
BB oscillations by plotting the so-called mixed fraction 

V, i 1 

FmiJt) = - ^— = - (1- cos Am*), (61) 

W Vos+Vls 2 

which exhibits a sinusoidal behavior with J r m i x (0) = 0. 

In summary, we would like to point out that a measurement of the time depen- 
dence of BB oscillations is similar to a B lifetime measurement, the difference being 
that the entire sample is divided into a mixed sample and an unmixed sample. In 
the following section we review some of the B°B° mixing results at CDF. These 
measurements are an excellent proof that B flavour tagging works in a hadron col- 
lider environment. Besides a measurement of Am they allow, at the same time, the 
determination of the effective tagging efficiency eD 2 of the applied flavour tag. The 
frequency of the oscillation will determine Am, while the amplitude characterizes T>. 



7.2. B°B° mixing using same side tagging 

The measmement of the B°B° oscillation frequency Am^, using a same side tagging 
techniqueE!rE2l (SST) is based on partially reconstructed semileptonic B decays to 
D^'ivX. The lepton charge tags the B flavour at decay time^j/hile a same side 
tag provides the B flavour at production. It has been suggested^ that the electric 
charge of particles produced near a B meson can be used to determine its initial 
flavour. This can be understood in a simplified picture of fragmentation as shown 
in Figure |37|a). For example, if a b quark combines with a u quark to form a 
B~ meson, the remaining u quark may combine with a d quark to form a 7r + . 
Similarly, if a 6 quark hadronizes to form a J3° meson, the associated pion would be 
a 7r~ . Another source of correlated pions are decays of the orbitally excited (L = 1) 
B mesons B**° -> B^+n- or B**+ -> B^°tt + (see also Sec. 7.2.1). Here, 

no attempt is made to distinguish the hadronization pions from those originating 
from B** decays. 
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Fig. 38. Left: Invariant mass distributions of the five B — > ' uX signatures used in the 
analysis: (a) D° -> K+tt' in D°l+X, (b) D~ -> K+tt-tt- in D~l+X, (c) D*~ -> D°w~ with 
D° — > i^+7r _ and D° — > i<" + 7r - 7r + 7r~ in D*~^ + X, (d) X + 7r — 7r^ minus A"+7r _ mass difference 
for D*~ — » D°ir* with D° — > A" + 7r — 7r° in l + X. (e) Measured asymmetries as a function 

of proper decay length for the decay signatures: (top), D~£ + (middle), and the sum of all 

three D*~ £ + (bottom). The dashed lines are the results of the fit. 



As in the B lifetime analysis with partially reconstructed B mesons (see Sec- 
tion 6.4.2), B candidates are reconstructed using the decay chains B° — ► 
with D~ —>■ K + 7t~tt~, or D*~ — ► D°ir~ followed by D° decaying to K + ir~ , 
K + ir~ tt + it~ , or K + tt^t: 1 where ir~ denotes the low-momentum (soft) pion from 
the D*~ decay. Charged B mesons are reconstructed through B + — > D a t + v, with 
D Q — ► K + ir~ , where the D° is required not to form a D*~ candidate with any 
other 7r candidate in the event. The mass distributions of the four decay signatures 
with fully reconstructed D mesons are shown in Fig. 38 i), b) and c), while the 
distribution of m(K / K / K> f ) — m(K / n) for the signature with D*~ — ► D n~, followed 
by D° — ► K + tt~tt° (the 7r° is not reconstructed) is displayed in Fig. [38|d). 

To select the SST pion, all tracks within a cone of radius 0.7 in Tytp-space, cen- 
tered around the direction of the B meson, approximated by p(£) +p(D) are consid- 
ered (see Fig. |37jb). SST candidate tracks should originate from the B production 
point (the primary event vertex), and are therefore required to satisfy do/ado < 3, 
where Ud is_the uncertainty on the track impact parameter do. String fragmenta- 
tion modelscJ indicate that particles produced in the b quark hadronization chain 
have low momenta transverse to the direction of the b quark momentum. We thus 
select as the tag the track that has the minimum component of momentum, p£ , 
orthogonal to the momentum sum of the track and the B meson (see Fig. |37|b). 
The tagging efficiency, e, is defined as the fraction of B candidates with at least one 
track satisfying the above requirements. It is measured as e ~ 70% independent of 
the decay signature used. On average, there are about 2.2 SST candidate tracks 
per B candidate. 
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For each of the five decay signatures, the B candidates are subdivided into 
six bins in proper decay length, ct, and the B-n combinations are classified as 
right-sign (RS: B + ir~ and B°tt + ) or wrong-sign (WS: B + ir + and B°tt~). In 
addition, the asymmetry in the RS and WS combinations is formed, A(ct) = 
{NRs(cb) — Nws(ct))/(NRs(ct) + Nws(ct))' For B + mesons, we expect an asymme- 
try independent of ct: A+(ct) — const. = V + . The dilution T> + is a direct measure 
of the SST purity, where (l + X>)/2 is the fraction of correctly tagged events. Due to 
B°B° mixing, Ao(ct) for the neutral B mesons will vary as a function of ct. It follows 
from Eq. ( |57|) that the asymmetry is expected to oscillate as Ao(ct) = Vq-cos Arrid t. 
Mistags result in a decrease of the oscillation amplitude by the dilution factor T>q. 
The asymmetry is measured as a function of the proper decay length ct for both 
B + and B° mesons, and fit with the expected time dependence, obtaining Am^, 
T>o, and T> + . 

The measured asymmetries A meas (ct) are shown in Fig. ^8^). If the £ + D° and 
£+£)(*)- signatures are pure signals of B + and B° decays, Amj could simply be 
extracted using the time-dependence of Ao(ct). However, the signatures are mix- 
tures of B + and B° decays, and thus A meas (ct) is a linear combination of the true 
asymmetries Ao(ct) and A+(ct). To extract Arrid, T>q and T> + , it is necessary to 
determine the sample composition of each D^£ + signature, which is the fraction of 
D^£ + candidates originating from the decays of the B° and B + mesons. Because 
a B + is associated with a 7r~, whereas an unmixed B° is associated with a tt + , the 
observed asymmetries are reduced by cross-contamination which can arise if the soft 
pion 7r~ from the D*~ decay is not identified. The decay sequence B° — ► D*~£ + v 
will be reconstructed as £ + D°, that is, as a B + candidate. Another source of 
cross-contamination arises from semileptonic B decays involving P-wave D** reso- 
nances as well as non-resonant D^tt pairs, which cannot be easily recognized and 
removed from the sample. For example, the decay sequence B° — > D**~ £ + v, fol- 
lowed by D**~ — > D°tt~^ (by 7r** we denote the pion originating from a D** decay) 
will be reconstructed as £ + D°, because of the missed 7r~„. Again, a B° decay is 
misclassified as a B + candidate. The tagging is further complicated when a tt^ 
from a £>** decay is present. The irf^ may be incorrectly selected as the SST pion, 
always resulting in a RS correlation. The requirement do/<7d < 3, described above, 
reduces this effect. The 7r„* originates from the B meson decay point, whereas the 
appropriate tagging track comes from the B meson production point. 

Taking these effects properly into account, the mass difference Ato^ and the 
dilutions T>q and T>+ are determined from a x 2 -ftt to the measured asymmetries 
A mcas {ct) with the fit result overlaid, as shown in Fig. ^8|s). The oscillation in the 
neutral B signatures is clearly present. The final result for the mixing frequency is 

Am d = (0.471 t° f 6 l ± 0.034) ps" 1 . (62) 

In addition, the following values for the neutral and charged SST flavour tagging 
dilutions are obtained: 

B° : X> = (18 ± 3 ± 2)%, eV\ = (2.4 ± 0.7l^)%, (63) 
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B+ : V+ = (27 ± 3 ± 2)%, eV\ = (5.2 ± 1.2±g;|)%. (64) 

The fit indicates that ~ 82% of the L>°^ + signature comes from B + decays, while 
~ 80% of the D~£+ and - 95% of the D*~£+ originate from B°. The B° component 
of the D°£ + signature can be seen as a small anti-oscillation in Fig. |38|e), top. 

7.2.1. Production of B** mesons 

To shed some light on the question of how much of the observed tagging power of 
same side tagging arises from fragmentation tracks compared to pions from B** 
decays, we briefly discuss the production of B** mesons at CDF. The term B** 
is a collective name for the four lowest lying L = 1 states of B mesons which are 
usually labeled by the total angular momentum J of the light quark resulting in two 
doublets with J = 1/2 and 3/2. The states in the J = 1/2 doublet are expected to 
be broad (a ~ 100 MeV/c 2 ) since they can decay through a S- wave transition while 
the J = 3/2 states decay through a D-wave transition and are therefore exnected 
to be narrow (a ~ 20 MeV/c 2 ). B** states have been observed at LEPo with 
properties in reasonable agreement with the expectations. 

The procedure for selecting B mesons follows the one described above in Sec. 7.2 
reconstructing semileptonic B decays into D^*'ivX. About 5500 D* + £ and D + £ 
candidates associated with a B° signature and about 4200 D°£ events indicating a 
B + signature are selected. B** candidates are constructed by combining the B can- 
didates with all tracks compatible with those originating from the primary interac- 
tion vertex. Since the resolution on the B** invariant mass is impaired by the un- 
known momentum of the missing neutrino, a quantity Q = m(Bir) — m(B) — m(w) 
is constructed. A resolution of ~ 50 MeV/c 2 on Q is expected from B** decays. 
Since the charge of the pion from a B** decay always matches the B flavour, the 
sample is divided into right-sign and wrong-sign combinations. We expect an excess 
of RS combinations confined to Q < 500 MeV/c 2 corresponding to the B** states. 

The essential task of this analysis it to properly determine the backgrounds to 
the B** signal. These can be divided into correlated and uncorrelated components. 
Three sources of uncorrelated background are taken into account: Combinatorial 
background from fake B meson candidates, particles from the underlying event 
and tracks from events with multiple hard pp collisions. All three components are 
measured from data and subtracted from the Q distributions. The main correlated 
background arises from fragmentation particles which contribute to the right-sign 
excess. The shape of this background is predicted using Monte Carlo calculations 
tuned to CDF data (see Sec. 4.5.2). To extract the B** production fraction from 
the Q distribution, cross contamination between B° and B + via D** and B* states 
is taken into account as well as the known flavour mixing of B a mesons. 

The result of the fit to the Q distributions is shown in Fig. |39|. The points are 
the data, the dashed curves are the fitted shapes of the hadronization component, 
the dotted histograms include all backgrounds, and the solid histograms are the 
sum including the fitted B** signal. The B** production fraction, defined as the 
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Fig. 39. Result of the fit to the Q distributions of the Bit candidates. The dashed curves are 
the fitted hadronization components, the dotted histograms include all backgrounds, and the solid 
histograms are the sum including the fitted B** signal. 



probability that a B quark hadronizes into a B** state, is measured as /** = 
0.28 ± 0.06 ± 0.03. The experimental resolution does not allow disentangling the 
four B** states and the average mass of the ensemble is varied collectively to obtain 
m(.Bi) = (5.71 ± 0.02) GeV/c 2 for m{B 1 ) of the narrowest state (J = 3/2). 

7.3. B°B° mixing using jet-charge and lepton tagging 

The analysis reported next@ uses the same single lepton data samples as the SST 
mixing result but increases the number of B mesons by over an order of magnitude 
by inclusively reconstructing B hadrons decaying semileptonically rather than using 
lepton-charm correlations. The inclusive reconstruction is based on identifying the 
B hadron decay point by associating the trigger lepton with other B decay products 
to form a secondary vertex. 

First, charged particle jets are reconstructed in the event using a track based jet 
clustering algorithm (see Sec. 4.4). The trigger lepton is associated with a jet. The 
search for the B decay point in the trigger lepto n je t is based on the technique to 
identify b quark jets coming from top quark decaya£3. Some modifications are made 
to maintain good efficiency for reconstructing displaced vertices, since B hadrons 
in this data sample have substantially lower p t than B hadrons from top decays. 
Tracks in the jet are selected for reconstructing the secondary vertex based on the 
significance of their impact parameter <io/0d o with respect to the primary vertex, 
where Ud a is the estimate of the error on do. First, displaced vertices containing 
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three or more tracks satisfying a loose set of track quality requirements are searched 
for. If no such vertices are found, two-track vertices which satisfy more stringent 
quality requirements are accepted and the two-dimensional decay length of the 
secondary vertex is calculated. To reduce background from false vertices formed 
from random combinations of tracks coming from the primary interaction vertex, 
|-L xy /(7L xy | > 2.0 is required, where <Jh xy is the estimated error on L xy . This leaves 
243,800 events: 114,665 from the electron data sample and 129,135 from the muon 
data sample. Secondary vertices are also searched for in the other jet in the event. 
If an additional displaced vertex is found, the event is classified as a "double- vertex" 
event. 

The flavour of the B hadron producing the trigger lepton is identified using a 
jet-charge tag and an opposite side lepton tag. First, an additional lepton (e or n) 
from the semileptonic decay of the opposite side B hadron is searched for in the 
event. The invariant mass of this lepton and the trigger lepton must be greater 
than 5 GeV/c 2 to reject leptons that come from the same B hadron producing the 
trigger lepton. Approximately 5.2% of the 243,800 events contain an opposite side 
lepton. If such a lepton is not found, the jet produced by the opposite b quark is 
identified by calculating a quantity called the jet-charge Qj C t, 

Q JCt = ( 6 5) 

where qi and pi are the charge and momentum of track i in the jet and a is a unit- 
vector defining the jet axis. On average, the sign of the jet-charge is the same as the 
sign of the b quark that produced the jet. In the case of a double- vertex event the jet 
containing the other secondary vertex is used to calculate Qj e t- About 7.5% of the 
sample consists of jet-charge double-vertex events and 42% are jet-charge single- 
vertex events. Figure [l0| shows, on the left side, the jet-charge distributions for 
single-vertex and double-vertex events. The degree of separation between the solid 
and dashed distributions is related to the tagging power of the jet-charge flavour 
tag. For double- vertex events, the presence of the second displaced vertex increases 
the probability that the selected jet in fact originates from the other B hadron in 
the event. This translates into a significantly larger separation in jet-charge for 
double-vertex events. 

Unlike the case of the SST mixing analysis, where a charm peak in the invari- 
ant mass spectrum indicates the amount of signal and background, the number of 
trigger leptons from bb production must be determined differently in this analysis. 
Other sources of the selected events include cc production and light quark or gluon 
production that can result in fake leptons or vertices. The fraction of events from 
false vertices and fake leptons is found to be small. The fraction of events due to 
bb and cc production is determined using two kinematic quantities: the trigger lep- 
ton p[ cl and the invariant mass m cl of the cluster of secondary vertex tracks. The 
quantity p[ el is defined as the magnitude of the component of the trigger lepton 
momentum perpendicular to the jet axis. The trigger lepton is removed from the 
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Fig. 40. Left: Jet-charge distributions for double-vertex (top) and single-vertex (bottom) events 
for the e- and /^-trigger data. Right: Same-sign fraction as a function of proper decay length. A 
representation of the fit result is superimposed on the data. 



jet, and the jet-axis is recalculated to determine p r t . To calculate m , the pion 
mass is assigned to all of the tracks used to form the secondary vertex (except the 
lepton) . The trigger lepton is included in this calculation even if it is not attached 
to the secondary vertex. Both kinematic quantities are effective in discriminating 
between bb and cc events because of the significant mass difference between hadrons 
containing b and c quarks. Template p\ cl and m cl distributions are obtained from 
bb and cc Monte Carlo samples (see Sec. 4.5.2) and fit to the data. The determi- 
nation of the sample composition yields the data being primarily (> 90%) from 
bb production. 

According to the flavour tag and the charge of the trigger lepton, the data 
are divided in a like-sign and opposite-sign sample. The proper decay time of the 
B hadron is determined from the L xy of the secondary vertex, combined with an 
estimate of the B momentum obtained from Monte Carlo calculations. From an 
unbinned maximum log-likelihood fit to the time dependence of the OS and LS 
sample, the B° oscillation frequency is determined to be 

Am d = (0.500 ± 0.052 ± 0.043) ps -1 . (66) 

The distribution of the like-sign fraction as a function of proper decay length is 
shown in Fig. [lO] on the right hand side. A representation of the fit result is super- 
imposed on the data. In addition, the effective tagging efficiencies of the jet-charge 
and opposite side lepton tags are obtained: 

jet-charge tag : eV 2 = (0.78 ± O.f 2 ± 0.08)%, (67) 

lepton tag: eV 2 = (0.91 ± O.f ± 0.11)%. (68) 

These values for eT> 2 are much lower than those achieved in experiments on the 
Z° resonances. On the other hand, the much higher bb cross section at the Tevatron 
can be used to compensate for the disadvantage in eT> 2 . 
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7.4. B°B° mixing using dilepton data 

The previous analysis searches single lepton events for additional leptons used as 
opposite side flavour tags. However, there exist large dilepton data sets at CDF 
where the lepton tag has already been triggered on. These data are also used for 
proper time dependent measurements of BB oscillations. In these data sets both 
leptons are assumed to come from the semileptonic decay of both B hadrons in the 
event: b — > t-yX and b — > I2X. This means, the flavour of the B meson at decay is 
tagged by its semileptonic decay, while the lepton from the semileptonic decay of 
the other B hadron in the event tags the B flavour at production. The requirement 
milxl-z) > 5 GeV/c 2 ensures that both leptons originate from two B hadrons and 
not from a sequential decay of one B hadron: b — > cl\X, with c — > l^X. Here, 
we briefly describe two mixing results using the e/z trigger data as well as dimuon 
trigger events. 

The analysis using e, dat E also searches for an inclusive secondary vertex 
associated with one of the leptons. As in the previous analysis in Sec. 7.3, secondary 
vertices are reconstructed using a modified version of the algorithm used to find 
displaced vertices in the search for the top quarkEJ. The decay length L xy of this 
vertex and the momenta of the tracks associated with the lepton provide an estimate 
of the ct of the B meson. Again, the important task of this analysis is to determine 
the sample composition, the fraction of events coming from bb decays with respect 
to events from cc or background events. The sample composition is estimated from 
several kinematic quantities, like p T t el or the invariant mass of the tagged secondary 
vertex. Here, pj; cl is defined as the transverse momentum of the lepton with respect 
to the highest pt track in a cone around the lepton. The determination of the sample 
composition finds that more than 80% of the events originate from bb decays. 

From a fit to the like-sign lepton fraction, as a function of ct, the mixing fre- 
quency Amj is extracted as shown in Fig. ffl|a). The fit includes components for 
direct and sequential b decays, cc, and fake events. In about 16% of the events with 
a secondary vertex around one lepton, a secondary vertex is also found around the 
other lepton. These events enter the like-sign fraction distribution twice, where we 
allow for a statistical correlation between the two entries. The final fit result is 

Am d = (0.450 ± 0.045 ± 0.051) ps -1 , (69) 

where the dominant systematic error arises from the uncertainty in the sample 
composition. 

The measurement of the B a B° oscillation frequency using dimuon dataE3 is very 
similar to the e/i analysis. Here, the algorithm used to find displaced vertices around 
one of the muons is based on the correlation between the impact parameter do and 
the azimuthal angle tp of tracks coming from these vertices. A vertex is formed 
by three or more correlated tracks with a significance do/(Jd > 2. Tracks from a 
secondary vertex form a line in the tp do-plane with non-zero slope, while tracks from 
the primary interaction vertex have small do values and show no correlation with (p. 
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Fig. 41. Fraction of (a) like-sign e/i events and (b) like-sign dimuon events as a function of 
pseudo-proper decay length. The results of the Araj fits are superimposed. 



The determination of the sample composition again uses different kinematic 
variables and finds the fraction of events from bb production to be larger than 80%. 
The final sample amounts to 2044 like-sign and 3924 opposite-sign dimuon events. 
The fraction of events with like-sign muons is shown in Fig. f4l|b) as a function of 
proper time of the B hadron associated with the identified secondary vertex. From 
a x 2 " n t to the like-sign fraction, the B° oscillation frequency is determined to be 

Am d = (0.503 ±0.064 ±0.071) ps"\ (70) 

where the main systematic error arises from the determination of the background 
and the fraction of muons from sequential charm decays. 

7.5. B°B° mixing using charm decays with a lepton tag 

We briefly describe two more B°B° mixing resultsll also using a lepton tag to iden- 
tify the B flavour at production. Both preliminary analyses allow the determination 
of the dilution or mistag probability of the lepton tag. One analysis uses higher 
momentum leptons than the other which can shed some light on the momentum 
dependence of the tagging dilution. 

The first analysis uses the dimuon as well as e/z trigger data and searches for 
a charm decay around one of the leptons identifying semileptonic B decays. In 
particular, leptons associated with a D* + meson are used as signature for B° decays. 
The D* + mesons are reconstructed in the decay mode D* + — > D°tt + , followed by 
D° — > K~ir + , K~ir + Tr + ir~ or K^tt + tt°. A signal of about 500 such decays is 
found on low backgrounds. The decay vertices of the D* + £~ combinations are 
reconstructed and the proper decay length of the B° meson is estimated using the 
D* + i~ momentum. The charge of the final state identifies the B° flavour at decay 
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Fig. 42. Asymmetry distribution of (a) D*~*~£~ candidates, (b) D~ candidates, and (c) D*~ 
candidates as a function of proper decay length. The points are the data and the curves show the 
result of the Am^ fits. 



time, while the other lepton in the event provides the B flavour at production. 
As in the B lifetime analysis using partially reconstructed B decays (Sec. 6.4.2), 
D* + l~ combinations are a good signature for neutral B mesons, but D** decays 
allow a contamination of charged B decays into the D* + l~ sample, diluting the 
charge correlation between the final states and the parent B meson. Using Monte 
Carlo techniques as in Sec. 6.4.2, the cross contamination from B + decays is found 
to be 0.19^;^ for the D*+ pT sample and 0.14±°;^ for the D* + e~ sample. 

The signal sample is divided in OS and LS events and the corresponding asym- 
metry distribution is shown in Fig. |42]a). From the time dependence of the OS 
and LS sample, the B° oscillation frequency Amj as well as the mistag probability 
Pw of the lepton tag are simultaneously extracted. The unbinned log-likelihood fit 
yields 

Am d = (0.512 ™ t°o° Q ll) ps" 1 , (71) 

while the flavour misidentihcation probability is pw = 0.302 ± 0.037^qq"27 yielding 
a dilution of the lepton tag of V = 0.396 ± 0.074lg$$. 

The second analysis uses the trigger lepton from the single lepton trigger data as 
B flavour tag. The transverse momentum of these leptons is greater than 8 GeV/c, 
in comparison to the dilepton data where pt(£) is greater than about 2 GeV/c. The 
B° meson is inferred from a decay opposite the trigger lepton, where the 

charge of the D^*' meson tags the B° flavour at decay. This means, this analysis 
searches for events with b — > (r X and b — > B° — > D^~X. A £~D^~ combination 
is therefore a signature for an unmixed event, while a £ + D^*' pair signals a mixed 
event. candidates are reconstructed as D~ — > K + tt~-k~ and D*~ — > D°tt~ , 

with D Q — > K + n~ . The decay length of the charm meson provides an estimate of 
the B° decay time together with the measurement of pt(D^~) used to infer p t (B°). 
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Fig. 43. Compilation of (a) the CDF time dependent B°B° mixing measurements and (b) the 
Anij average results for the LEP, SLD, and CDF experiments as of July 1998. 



Signals of (460 ± 31) D~ events and (358 ± 32) D*~ events are selected. The 
mesons may originate from bb and cc production. In addition, cross contami- 
nation of Z)'*) - mesons from charged B mesons via D** decays must be considered, 
as discussed above and in Sec. 6.4.2. Using MC techniques and fits to kinematic 
variables, it is found that approximately 65% of the signal sample comes from 
B° decays, ~20% from B + , and ~ 15% from cc production. The time dependent 
asymmetry distribution of £ + D^~ and l~ pairs is plotted in Fig. |42)b) for 
D~ (top) and D*~ (bottom). The fit for the B a mixing frequency results in 

Am d = (0.562 ± 0.068 ±°:o5o) ps" 1 , (72) 

while the dilution of the lepton tag is determined to be V = 0.572 ± O.OSOt^osi- 

7.6. Summary of B°B° mixing measurements 

A compilation of all Amj measurements at CDF is found in Figure |43|a) . We use 
these measurements to determine a CDF average result for Amd- Since the mea- 
surements of Amd are quite precise, it is important to correctly take the correlated 
systematic uncertainties into account in forming the average. In addition, most of 
the presented results depend on common physics parameters, where different values 
are used in the original analyses. In order to correctly handle these issues, we use 
the procedure developed by the LEP B Oscillation Working Group and adjust the 
measurements of Amd on the basis of a common consistent set of input values. 
First, the procedure linearly rescales each measurement of Amd for each common 
physics parameter, in accordance with the difference between the originally used 
parameter value and its new common value, and with the corresponding system- 
atic uncertainty determined in the analysis. This systematic uncertainty is also 
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adjusted, if the uncertainty on the parameter is different from that originally used. 
The combination procedure then performs a common fit of Ato^ and of the physics 
parameters related to the sources of the common systematic uncertainties. We also 
take into account a statistical correlation of about 20% between the Am^ result 
using jet-charge and lepton tagging and the Amj analyses using dilepton data. We 
determine the average of all Amj results at CDF to be 

Am d = (0.494 ± 0.026 ± 0.026) ps" 1 . (73) 

A comparison of this CDF average with the Amd averages of the LEP and SLD 
experiments as of July 1998 is displayed in Figure [43|b). The CDF Am^ results are 
quite competitive with other experiments. 

8. Search for CP Violation 

In 1964, Christenson, Cronin, Fitch, and Turlay discovered^ that the laws of nature 
are not invariant under the combined action of charge conjugation C and parity P. 
They found that the K\ meson, associated with the CP-odd state in the neutral 
kaon system, decays with a branching fraction of about 2 • 10~ 3 into 7r + 7r~, which 
is a CP-even state. This result can briefly be illustrated in the following way: As 
in the B meson system, discussed in the Sec. 6.5.3 and Sec. 7.1, neutral kaons exist 
in two CP-conjugate flavour states \K°) = \sd) and \K°) = \sd). Applying the CP 
operation to these states, we define the following phase convention: 

CP\K°) = -\K°) and CP \K°) = —\K°). (74) 

The physical kaon states with well defined mass and lifetime, which decay through 
the weak interaction, are linear combinations of the two CP-conjugate quark states: 

\K° 1 ) = 1/V2(\K°)~\K )) and \K° 2 ) - l/y/2 (\K°) + \K )). (75) 

If CP is conserved, the physical kaon states are CP eigenstates: 

CP\K°) = CP [1/V2(|A'°) - \K ))} = 1/V2{-\K°) + \K )) = +\K°). (76) 
We call Kl the "CP-even" state. For K$, we obtain 

CP\K° 2 ) = CP[1/V2(\K°) + \K ))} = 1/V2(-\K°) - \K )) = -\K°) (77) 

and call K® the "CP-odd" state. From experiment, we know that neutral kaons 
decay either into two or three pions. Since CP \it + it~) = +1 and CP |7r + 7r~7r°) = 
-1, we associate \Kf) with K% -> tt+tt" and with K° L -> tt+tt~tt°. If CP is 
conserved, a K^ L is not allowed to decay into two pions. The discovery of K\ decays 
into 7r + 7r~ thus manifests the existence of CP violation and states that the laws of 
nature are not exactly the same for matter and antimatter. 

More than 30 years have passed since the surprising discovery of CP violation, 
and the kaon system is still the only place where CP violation has been observed 
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in nature. However, the system of neutral B mesons is expected to yield large 
CP violating effects. This is the subject of this section, where we report on a 
search for CP violation in B° — ► J/tpKg decays at CDF. First, we briefly introduce 
the Cabibbo-Kobayashi-Maskawa mixing matrix and the CKM unitarity triangle. 

8.1. CKM matrix and unitarity triangle 

CP violation in the kaon system can be explained in an elegant way in the Standard 
Models with three generations, as originally suggested by Kobayashi and Maskawafl 
in 1973. Note, this was still before the discovery of the J/ip resonanceEj in 1974. In 
the Standard Model with SU(2)xU(l) as the gauge group of electroweak interac- 
tions, the quark mass eigenstates are not the same as the weak flavour eigenstates. 
The matrix which relates the mass eigenstates to the weak eigenstates was defined 
for six quarks and given an explicit parametrization by Kobayashi and Maskawacl. 
It generalizes the four-quark case, where the corresponding matrix is parametrized 
by a single parameter, the Cabibbo anglefl. The mixing is often expressed in terms 
of a 3 x 3 unitarity matrix Vckm , called the Cabibbo-Kobayashi-Maskawa mixing 
matrix, which operates by convention on the charge —1/3 quarks d, s, and b: 

d! \ ( V ud V us V ub \ / d \ 

s' = V cd V cs V cb s . (78) 

V J V V u V u Va, J \b J 

The individual matrix elements can, in principle, all be determined from weak 
decays of the relevant quarks or from deep inelastic neutrino scattering. 

There are several parametrizations of the CKM matrix. In the standard para- 
metrization, three rotation angles and an imaginary phase must be used. This phase 
being non-zero gives rise to CP violation in the weak interaction. Another popular 
approximation of the CKM matrix which emphasizes the hierarchy in the size of 
the matrix elements is, due to WolfensteinEa 

/ l-A 2 /2 A AX 3 (p-iri) \ 

Vckm=\ -A l-A 2 /2 AX 2 +G(A 4 ). (79) 

\ A\ 3 (l - p- irf) -AX 2 1 / 

Here, the CKM matrix is parametrized in powers of A, the sine of the Cabibbo 
angle A = sin# c = \V US \. The Cabibbo angle is well measured from leptonic kaon 
decays K — > -kIv yielding A = 0.2196 ± 0.0023. The other parameters used in the 
Wolfenstein parametrization are A, p, and rj. The parameter A is related to the 
CKM matrix element V c b = AX 2 and can be determined from semileptonic b — > c 
transitions. As we will see below, the parameters p and 77 are related to CP violation 
in the B meson system and are connected to the matrix elements V u b and Vtd- 

Because of the unitarity of the CKM matrix derived from conserving probabili- 
ties, we obtain a relation involving the two smallest elements V u b and Vtd of Vckm 



Vud.v* b + v cd v; b + V td V t * b - 0. 



(80) 
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Fig. 44. The CKM unitarity triangle (a) in standard representation and (b) in the (p, r;)-plane 
using the Wolfenstein parametrization. 



This is a triangle relation in the complex plane as shown in Fig. |44|a). The Cabibbo- 
Kobayashi-Maskawa matrix elements determine the size of the legs of this triangle, 
often called the CKM unitarity triangle. The angles a, ft, and 7 are related to 
CP violating asymmetries in B decays. For example, an asymmetry in B a — > 
J/ipKg decays measures sin 2ft, while B° — > tttt decays are related to sin 2a. For 
CP violation to be permitted in the Standard Model, the area of the unitarity 
triangle must be non-zero and the angles of the triangle different from zero or n. 

We briefly illustrate the current knowledge about the CKM unitarity triangle. 
For this purpose, we use the Wolfenstein parametrization, where the unitarity tri- 
angle is described in the (p, ?7)-plane, as shown in Fig. |44|b). We make the following 
approximations \V t * b \ = 1, V c d = —A, and set V u d = 1 — A 2 /2 1. We also normalize 
the legs of the triangle by \V c b\ keeping the base of the triangle of unit length. The 
other two legs are then described by \V* b /XV* b \ and \V t d /AV c * b |. A measurement of 
IKtft/Vcfel constrains one leg of the triangle to lie between the two half circles centered 
at the origin, as shown in Figure ^5|. The other leg of the triangle is constrained by 
our current knowledge of Vtd obtained from B°B° mixing measurements. The large 
uncertainty with which the B meson decay constant fg and the bag parameter (see 
Sec. 7.1) are known, allows this leg of the unitarity triangle to lie in the broad area 
between the two circles centered around the point (1,0). Finally, the knowledge 
of CP violation in the kaon system adds another constraint on the CKM unitarity 



triangle indicated by the third band in Fig. 45. We fit the three bands to originate 
from a common area and obtain the 95% confidence level contour, also shown in 
Fig. [|5| We see that the apex of the unitarity triangle, described by the coordinates 
(p, ?y), can be placed in a fairly wide range. Direct measurements of the angles a, 
ft, or 7 would better constrain the CKM unitarity triangle. 

8.2. Measurement of CP violation parameter sin 2/3 in B° — > J/xjjK^ 

One way to observe CP violation in B meson decays is to use the interference 
between the direct decay of a B° into a CP final state fcp and the process B° — > 
B° — > fcp which can provide a second amplitude to interfere with the direct decay 
B° — > fcp- A large CP violating effect is expectedLJ in B°/B° decays to the 
CP eigenstate J/tpKg. The interference of direct decays, B° — > J/ifjKg, versus 
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those that undergo mixing, B° — > B° — > J/ipK s , gives rise to a decay asymmetry 

Acp(t) = = sin2/? ■ sin Amdi; (81) 

where B°(t) (B°(t)) is the number of decays to J/ipK s at proper time t given that 
the produced meson was a J3° (B°) at t = 0. The CP phase difference between 
the two decay paths appears via the factor sin 2/3, and the BB flavour oscillation 
through the mass difference Am^ between the two B° mass eigenstates. Figure [46|a) 
illustrates the connection between the time dependence of the B°B° mixing asym- 
metry Anix (top) and the CP asymmetry Acp (bottom) in e.g. B° — > J/ipK s 
decays. At t = we start with a pure B° state. Since no B° mesons are present to 
interfere with, the CP asymmetry is zero. After about two B lifetimes, we have the 
same amount of B° and B° due to mixing (.4 m i x = 0) and the CP asymmetry is 
maximal. From this we see that the CP asymmetry follows a sinusoidal time behav- 
ior with sin 2/3 characterizing the amplitude of this sine curve as shown in Eq. (pi]) 
while the mixing asymmetry follows a cosine like time curve as expressed in Eq. ( |57| ) . 

Themeasurement of the CP violation parameter sin 2/3 using B° /B° — > J/ipK s 
decaysES reconstructs J/ip mesons through and searches for Kg — > 7r + 7r _ 

candidates. The J/ip and Kg daughter tracks are combined in a four particle fit to 
originate from a common B° vertex. The decay length of the B° is used to calculate 
its proper decay length ct. We define Mn = (wfit — ma)/ & fit, where mat is the mass 
of the B candidate from the fit described above, crgt is its uncertainty (typically 
~ 9 MeV/c 2 ), and too is the central value of the B° mass peak. The normalized 
masses Mn for the selected candidates with ct > are shown in Fig. fltfb) along 
with the result of the maximum log- likelihood fit described later. The fit yields 
(198 ± 17) B°/B° mesons for all ct. 

To measure Acp{t), we need to know whether the production flavour of the 
B meson is B° or B°. We determine this with a same side tagging method described 
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Fig. 46. (a) Illustration of the connection between the time dependence of the B°B° mixing 
asymmetry (top) and the CP asymmetry (bottom) in B° —> J/tpK^ decays, (b) Normalized mass 
distribution of J/ipK^ candidates. The curve is the Gaussian signal plus linear background. 



in Sec. 7.2. The effectiveness of this method is demonstrated by tagging B — > D^tv 
decays and observing the time dependence of B°B° oscillations measuring Am^. 
Applying the SST method to the J/ipKg sample yields a tagging efficiency of ^65%. 
Since negative (positive) tags are associated with J5°'s (B°'s), the asymmetry 

a, \ N-(ct)-N+(ct) 



is formed analogous to Eq. (pi]), where N ± (ct) are the numbers of positive and 
negative tags in a given ct bin. The sideband-subtracted asymmetry is displayed in 



Fig. 47 1) along with a x 2_n t (dashed curve) to Aq sin Am^t, where Am^ is fixedH 
to 0.474 ps _1 . The amplitude, Ao — 0.36 ± 0.19, measures sin2/3 attenuated by 
the dilution factor T>q = 2pn — 1, where pu is the probability that the tag correctly 
identifies the B flavour. 

The fit is refined using an unbinned maximum log-likelihood method. This fit 
makes optimal use of the low statistics by fitting signal and background distributions 
in Mjq and ct, including sideband and ct < events which help constrain the 
background. The likelihood fit also incorporates resolution effects and corrections 
for systematic biases, such as the inherent charge asymmetry favouring positive 
tracks resulting from the wire plane orientation in the main drift chamber. The 
solid curve in Fig. [47|a) is the result of the likelihood fit, which yields 2?osin2/3 = 
0.31 ± 0.18± 0.03 including systematic effects. As expected, the two fits give similar 
results, indicating that the result is dominated by the sample size and that the 
corrections and improvements of the likelihood fit introduce no dramatic effects. 
Also shown in the Fig. f47|a) inset is the relative log-likelihood as a function of 
I\>sin2/3. The shape is parabolic, indicating Gaussian errors. 
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Fig. 47. (a) The sideband-subtracted tagging asymmetry as a function of the reconstructed 
J/tpK^ proper decay length (points). The dashed curve is the result of a simple x 2_n t to 
.ArjsinAm^t. The solid curve is the likelihood fit result, and the inset shows a scan through 
the log-likelihood function as T>q sin 2/3 is varied about the best fit value, (b) "Confidence belt" of 
the true value of sin 2/3 versus the measured value of sin 2/3 used to extract the limit on sin 2/3. 



To obtain sin 2/3, dilution measurements from other B samples^are used. CDF's 
best single V measurement from a large B -> D^IX sampleBE3 is 0.18lt°;°32 
(see Sec. 7.2). Because of differing lepton p t trigger thresholds, the average p t of the 
semileptonic B sample is ~21 GeV/c, but it is only ~ 12 GeV/c in the J/ipKg data. 
This difference is corrected for using Monte Carlo studies. The T>o appropriate for 
the J/ipKg sample is found to be 0.166 ± 0.018 ± 0.013, indicating a small shift 
from 0.181. The first error is due to the uncertainty in the dilution measurement, 
and the second is due to the Monte Carlo extrapolation. 

Using this value of T>q, sin 2(3 is determined to be 

sin 2/3 = 1.8 ±1.1 ±0.3. (83) 

The central value is unphysical since the amplitude of the measured asymmetry is 
larger than T>q. To express this result in terms oijl confidence interval in sin 2(3, the 
frequentist construction by Feldman and CousinsE3 is followed. This approach gives 
proper confidence intervals even for measurements in the unphysical region. Fig- 
ure E7b) shows the "confidence belt" of the true value of sin 2(3 versus the measured 
value of sin 2(3 using the frequentist method. This measurement thereby corre- 
sponds to excluding sin 2(3 < —0.20 at 95% confidence level. This result favours 
current Standard Model expectations of a positive value of sin 2/3 and establishes 
the feasibility of measuring CP asymmetries in B meson decays at a hadron collider. 

8.2.1. Updated measurement of CP violation parameter sin 2(3 in B° — > J/ipKg 

Just prior to completion of this article, the CDF collaboration released an updated 
measurement of the CP violation parameter sin 2(3 in B° — > J/tpKg decays. This 
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Fig. 48. (a) Normalized mass distribution of J/ipKa candidates. The curve is the Gaussian signal 
plus linear background from the log-likelihood fit. (b) The true asymmetry sin 2/3 as a function of 
the reconstructed J/ipK^ proper decay length. The data points are sideband subtracted and are 
combined according to the effective dilution for single and double tags. The non-SVX events are 
shown on the right. 



preliminary result combines the same side tagging method with a lepton and jet 
charge flavour tag (see Sec. 7.3). In addition, it includes J/ipKg events which are 
not fully contained within the acceptance of the SVX detector. 

The selection criteria for the B° — > J/tpKg sample with J/tp — > and 
Kg — > 7r + 7r~ are very similar to the requirements used in the analysis in Sec. 8.2. 
To increase the sample size, both muon candidate tracks are no longer required to be 
measured in the silicon vertex detector. The data are divided into two samples. The 
non-SVX sample contains events where both muons are not required to have SVX 
information, accepting events that do not have a precise decay length measurement. 
However, about 30% of the events in this sample have one muon track with SVX 
information. The other sample of J/ipKg decays, called SVX sample, requires both 
muons to be well measured in the SVX providing precise decay length information. 
This sample is very similar to the data used in the sin 2(3 analysis described in 
Sec. 8.2. The normalized mass distribution (see Sec. 8.2) for all J/tpKg candidates 
is shown in Fig. [48|a) comprising a sample of (395 ± 31) signal events. The SVX 
sample contains (202 ± 18) signal events while the non-SVX sample consists of 
(193 ± 26) events. 

This analysis combines all three flavour tagging methods studied at CDF. In 
addition to the same side tagging algorithm, the two opposite side tagging methods, 
jet charge tagging and lepton tagging, as described in Sec. 7.3, are added. The 
tagging dilutions and efficiencies for both opposite side tags are determined from a 
sample of (985 ± 46) B+ -> J/^K+ events and are presented in Table @. The SST 
tagging dilution for the SVX sample (D = (16.6 ± 2.2)%) is taken from the sin 2(3 
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Table 7. Summary of efficiency e and dilution T> of the tagging algorithms. The quoted efficiencies 
are relative to the entire J/ipK^ sample. To compare the SST efficiency with the number given 
in Sec. 8.2, it is necessary to double the SST efficiency for the SVX sample. 



Flavour tag 


Data 


Efficiency e 


Dilution T> 


eV 2 


Same side tag 

Lcpton tag 
Jet charge tag 


SVX sample 
non-SVX sample 
all events 
all events 


(35.5 ± 3.7)% 
(38.1 ± 3.9)% 
(5.6 ± 1.8)% 
(40.2 ± 3.9)% 


(16.6 ±2.2)% 
(17.4 ±3.6)% 
(62.5 ± 14.6)% 
(23.5 ±6.9)% 


(1.0 ±0.3)% 
(1.2 ±0.5)% 
(2.2 ± 1.2)% 
(2.2 ± 1.3)% 



analysis summarized in Sec. 8.2. For the non-SVX sample, the SST algorithm is 
slightly modified to include non-SVX tracks as candidate tagging tracks. A dilution 
scale factor, which relates the SST tagging performance from the SVX sample to the 
non-SVX sample, is derived from the J/ifjK + data resulting in T> = (17 .4 ± 3.6)% 
for the non-SVX sample (see also Tab. f7|). 

Each event can be tagged by as many as two tags: One same side tag and one 
opposite side tag. If both the lepton and jet charge tags are available, only the 
lepton tag is used, following the B°B° mixing analysis described in Sec. 7.3. The 
procedure used to combine double tagged events calculates a combined dilution 
weighted by the individual dilutions and combines the efficiencies in a similar way. 
Defining the two tags qi and q 2 as +1 if they identify a B° meson, as —1 if they tag 
a B° and as if the tag is not applicable, the individual tags are weighted by the 
dilutions T>\ and T>2 according to T)' x = q\T)\ and T>' 2 = 52^2, respectively. With 
the individual efficiencies £\ and £2, the combined dilution T> com and efficiency e com 
are defined as 

^'1+^2 



V[V 2 



and s com = £i£ 2 (1 + T>'{D' 2 ), (84) 



where the sign of T> com is the combined tag with dilution |2? com |. 

In a similar way as summarized in Sec. 8.2, a maximum log- likelihood method 
is used to make optimal use of the statistics by fitting signal and background distri- 
butions in normalized mass and ct, including sideband events. The true asymmetry 
as a function of the reconstructed J/ipKg proper decay length is shown in Fig. |4g|b) 
separately for the SVX and non-SVX sample with the result of the fit overlaid. 
The non-SVX sample contribution is included as a single point because of the low 
decay length resolution. The full log-likelihood fit uses both the SVX and non-SVX 
sample and properly treats the decay length and error for each event. The final fit 
yields 

sin 2/3 = 0.79 ±0.39 ±0.16, (85) 

where the systematic error reflects the uncertainty in the dilution parameters. Al- 
though the individual dilutions are not precisely determined due to the limited 
statistics of the B + —* J/ipK + sample, this error does not dominate the uncertainty 
on sin 2/3. Removing the constraint that fixes Arrid to its world average valuetl, 
the fit determines sin 2/3 and Am ( i simultaneously to be sin 2/3 = 0.88lg'4^ and 
Amj = (0.68 ± 0.17) ps _1 as indicated by the dashed line in Fig. |^b). The result 
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given in Eq. ( |85| ) leads to the confidence intervaLsin 2/3 < —0.08 (95% C.L.) using 
the frequentist method by Feldman and CousinsE3. This result excludes sin 2/3 < 
at the 93% confidence level and is the best direct hint for CP violation in the neutral 
B meson system to date. 

9. Future B Physics at CDF 

The Fermilab accelerator complex is currently undergoing an upgrade to produce an 
order of magnitude higher luminosities in the Tevatron collider. The largest change 
is the replacement of the Main Ring with the new Main Injector, which will provide 
higher proton intensity onto the antiproton production target, and larger aperture 
for the antiproton transfer into the Tevatron increasing the antiproton current and 
thus the luminosity of the Tevatron. After the completion of the Main Injector, 
the Tevatron is scheduled to run again in April 2000 at a centre-of-mass energy of 
2.0 TeV. In this so-called Run II, luminosities of 2.0 • 10 32 cm _2 s _1 will be reached 
yielding an integrated luminosity of 2 fb^ 1 delivered to the collider experiments 
within two years. To handle the higher luminosities and shorter bunch spacing of 
initially 396 ns and later 132 ns, the CDF experiment is currently undergoing a 
major detector upgrade. In this section, we first briefly describe the CDF detector 
upgrade and then summarize the prospects for measuring sin 2/3 and sin 2a at CDF 
in Run II. 

9.1. CDF II detector upgrade 

The CDF detector improvements for Run II are motivated by the shorter accelerator 
bunch spacing and the increase in luminosity by an order of magnitude. The primary 
upgrade goal is to maintain detector occupancies at Run I levels, although many of 
the detector changes also provide qualitatively improved detector capabilities. The 
CDF II upgrade is described in detail elsewhereEll. A schematic view of the CDF II 
detector is shown in Figure f49]a). We briefly summarize here the changes most 
relevant for B physics in Run II. 

One major upgrade is to the charged particle tracking system (see Fig. |49)b) , vital 
for the CDF B physics program. A new silicon vertex detector (SVXII) will consist 
of five layers of double sided silicon from radii of 2.9 cm to 10 cm arranged in five 
axial layers, two small angle (1.2°) and three 90° stereo layers. SVXII will consist 
of three modules covering the entire pp luminous region and provide stand alone 
3-dimensional tracking. The silicon sensors are read out by a radiation hard chip 
in deadtimeless mode. In addition, an intermediate silicon layer (ISL) consisting of 
two double sided silicon layers at larger radii permits stand alone silicon tracking 
out to |?7| < 2. The inner layer provides complete rj coverage while the outer layer 
has partial coverage (1 < |?/| < 2). A new open cell drift chamber (COT) will 
operate at a beam crossing time of 132 ns without having overlapping events in a 
single cell. This is accomplished with a fast gas and shorter drift cells of 0.9 cm. 
The COT consists of 96 layers arranged in four axial and four stereo superlayers. It 
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Fig. 49. (a) Schematic cut-away view of the CDF II detector, (b) Longitudinal view of the CDF II 
tracking system. 



increases the number of stereo layers from 24 in the CTC to 48 and provides dE/dx 
information for particle identification. 

The upgrades to the muon system almost double the central muon coverage. 
The CMP coverage increases by ~ 17% while the CMX system is completed, in- 
creasing its coverage by ~45%. A new intermediate muon system (IMU) extends 
the muon coverage up to | rj \ < 1.5 with fine granularity and provides coverage suffi- 
cient to identify isolated high p t tracks as muons or hadrons between \r)\ of 1.5 and 
2.0. In addition, a new scintillating tile plug calorimeter will allow good electron 
identification up to \rj\ < 2. 

All front-end electronics is designed to handle the 132 ns beam crossing period. 
The trigger and data acquisition upgrade allows for higher data rates such as a 
50 kHz Level 1 accept rate and increases the sophistication of the trigger decision. 
Data are stored in a 42 cell pipeline while awaiting the Level 1 trigger decision and 
can be transfered to Level 2 with no deadtime to the Level 1 trigger. Finally, the 
CDF collaboration has recently proposed two "beyond the baseline" projects which 
will significantly enhance the B physics capabilities of the CDF II detector. These 
include the installation of a low-mass radiation hard single sided silicon detector, 
with axial strips at a very small radius of ~ 1.6 cm, just outside the beam pipe 
as well as the installation of a time-of-flight (TOF) system, employing 216 three 
meter long scintillator bars with fine mesh photomultiplier tubes on each end, to 
be located between the outer radius of the COT and the superconducting solenoid 
magnet. 

9.2. Measurement of sin 2/3 in Run II 

CDF has the advantage of being an existing experiment with plenty of data and 
experience from Run I. We will make use of the Run I knowledge to estimate the 
Run II expectations for measuring CP violation. For the measurement of sin 2/3 in 
the B° -> J/ipK° s channel, CDF expects 10,000 J/tpK a s events with J/tp -> n+a~ 
and Kg — > 7r + 7r~. This number is estimated in the following way: Starting with 
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Table 8. Summary of flavour tagging methods used in the measurement of sin 2/3 in Run II and 
the data samples used to calibrate the tagging algorithms. 



Flavour tag 


eV 2 


Calibration sample 


Sample size 


Same side tag 


2.0% 


J/i>K*° 


~ 20, 000 


Jet charge tag 


3.0% 


J/ipK+ 


~ 40, 000 


Lepton tag 


1.7% 


J/ipK+ 


~ 40, 000 



~ 200 J/tpKg events reconstructed in the SVX in Run I and multiplying it by 
2 fb _1 /110 pb _1 for the total Runll integrated luminosity, by a factor of 1.5 for 
the extended coverage of SVX II, and by two for a lower p t threshold of 1.5 GeV/c 
in the muon trigger, as well as an improved muon coverage with the completed 
CMX, one obtains 10,800 reconstructed J/tpKg events. CDF also plans to trigger 
on J/tp —> e + e~ which would increase the number of J /^K^ events by ~ 50%. 
However, these data are not included in this estimate. 

In Run II, CDF expects to improve the effective tagging efficiencies eV 2 of the 
B flavour tagging methods, as summarized in Table ||. The extended lepton coverage 
with the completed CMX, IMU, and the plug calorimeter results in a total eD 2 of 
1.7% for lepton tagging. A significant improvement in eT> 2 ~ 3% is possible for 
jet charge tagging. The extended coverage of SVX II and ISL as well as their 
added pattern recognition capabilities will substantially enhance the purity of the 
algorithm. Finally, for our Run II extrapolation a value of eV 2 ~ 2% is assumed for 
same side tagging. 

To estimate the uncertainty on a sin 2(3 measurement in Run II, we extrapolate 
from the measured error of the Run I result given in Eq. (|S3| ) which uses SVX events 
only. Defining the CP asymmetry as A = T> sin 2/3, we obtain sin 2/3 = A/T> and 
write the error on sin 2(3 as 

- 2 (sin2/ 3 )^(^) 2 + (sm2^) 2 , (86) 

where we substitute ^4/2? with sin 2(3 in the second term. This way, the error 
on sin 2/3 can be broken up into a statistical part from the measurement of the 
CP asymmetry A and a systematic part due to the dilution uncertainty. In Run II, 
the flavour tag dilutions will be calibrated with large samples of about 40,000 
J /ipK + and ~ 20, 000 J/i(K* Q events. Assuming sin 2/3 = 1 in the second term 
of Eq. ( |S6| ) as worst case and the same signal to noise ratios for the J/ipK samples 
as in Run I, the error on a measurement of sin 2/3 can be estimated to a 2 (sin 2/3) = 
0.078 2 + 0.031 2 . This will allow the observation of CP violation in Runll and a 
measurement of sin 2(3 from B° — > J/t(Kg with a precision of ±0.08 comparable to 
e + e _ machines. 

9.3. Measurement of sin 2a in Run II 

Another goal of B physics in Runll is the observation of CP violation in B° — > 
7r + 7r~ measuring sin 2a. The key to measuring the CP asymmetry in B° — > ir + ir~ 
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is to trigger on this hadronic decay mode in pp collisions. CDF plans to do this 
with its three level trigger system. On Level 1, two oppositely charged tracks with 
Pt > 2 GeV/c found with a fast track processor yield an accept rate of about 22 kHz 
at luminosities of ~ 1.0- 10 32 cm -2 s _1 . This rate will be reduced to less than ~ 25 Hz 
on Level 2 using track impact parameter information (do > 100 fmx). On Level 3, 
the full event information is available further reducing the trigger rate to about 
1 Hz. With this trigger, CDF expects about 10,000 B° -> tt+tt" events in 2 fb -1 , 
assuming BR(B° — > 7r + 7r~) = 1.0 • 1CP 5 . With the same effective tagging efficiency 
as in the sin 2(3 measurement above (see Tab. ||), CDF estimates an uncertainty of 
0.10 on sin 2a. Backgrounds from B — > Kir and B — > KK decays can be extracted 
from the untagged signal by making use of the invariant mass distribution, as well 
as CDF's dE/Ax capability with the COT. 

10. Conclusions 

The CDF collaboration has shown that it is possible to take advantage of the high 
production rate for B hadrons at the Tevatron Collider. The key elements for 
a broad B physics program at a hadron collider are the successful operation of 
a silicon vertex detector to identify B decay vertices displaced from the primary 
pp interaction vertex, the excellent tracking capabilities of CDF's central tracking 
chamber, together with the silicon vertex detector and the use of specialized triggers. 
Due to these features of the CDF detector together with the high yield of B hadrons, 
the hadron collider B physics program complements that at e + e~ machines. 

In this article, we gave a brief overview of heavy quark production in pp collisions 
and described some of the features of B physics at a hadron collider. We reviewed 
the B hadron lifetime measurements at CDF, all very competitive with results from 
LEP and SLC. We discussed several proper time dependent measurements of B°B° 
oscillations as well as B tagging methods to identify the B flavour in hadronic 
collisions. The application of all B flavour tags to the current data set of B° — > 
J/ijiKg decays at CDF established the feasibility of measuring CP asymmetries in 
B meson decays at a hadron collider. With many years of experience in B physics at 
the Tevatron, the CDF collaboration plans to observe CP violation in B° — ► J/ ipKg 
in the next run of the Tevatron collider scheduled to start in the year 2000 and will 
measure sin 2/3 with a precision comparable to the dedicated e + e _ B factories. 
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